chanter 0.0.9__tar.gz → 0.1.0__tar.gz

chanter-0.1.0/PKG-INFO

@@ -0,0 +1,25 @@
+Metadata-Version: 2.4
+Name: chanter
+Version: 0.1.0
+Summary: Simple quiescent galaxy spectral modelling and fitting
+Author: Struan Stevenson
+Author-email: struan.stevenson@ed.ac.uk
+Description-Content-Type: text/markdown
+Requires-Dist: numpy
+Requires-Dist: astropy
+Requires-Dist: matplotlib
+Requires-Dist: spectres
+Requires-Dist: nautilus-sampler
+Requires-Dist: json
+Dynamic: author
+Dynamic: author-email
+Dynamic: description
+Dynamic: description-content-type
+Dynamic: requires-dist
+Dynamic: summary
+
+# CHANTER
+
+#### CHANTER is a simple quiecent galaxy spectral energy distribution (SED) modelling/fitting code, based on the larger BAGPIPES code, of astronomical literary fame. CHANTER provides the barebones of Stellar Population Synthesis (SPS), as well as redshift and dust/IGM attenuation effects, and is written in a simple and digestible format, ideal for learning the theory behind larger spectral fitting software codes. CHANTER does not yet include nebular attenuation, and therefore only models the SED of quiescent galaxies.
+
+Currently, I cannot include the SSP/IGM grids in the package due to PyPi size limitations, but this is being appealed.

chanter-0.1.0/README.md

@@ -0,0 +1,5 @@
+# CHANTER
+
+#### CHANTER is a simple quiecent galaxy spectral energy distribution (SED) modelling/fitting code, based on the larger BAGPIPES code, of astronomical literary fame. CHANTER provides the barebones of Stellar Population Synthesis (SPS), as well as redshift and dust/IGM attenuation effects, and is written in a simple and digestible format, ideal for learning the theory behind larger spectral fitting software codes. CHANTER does not yet include nebular attenuation, and therefore only models the SED of quiescent galaxies.
+
+Currently, I cannot include the SSP/IGM grids in the package due to PyPi size limitations, but this is being appealed.

{chanter-0.0.9 → chanter-0.1.0}/chanter/galaxy_fitter.py

@@ -18,8 +18,10 @@ class galaxyfitter(object):
             Free parameters.
         """
         model = modelgalaxy(fit_instruc)
-        wavs, flux = model.get_spectrum()
-        effwavs, phot = model.get_photometry(wavs, flux, filt_files)
+        self.wavs, flux = model.get_spectrum()
+        self.spectrum_samples.append(flux)
+        effwavs, phot = model.get_photometry(self.wavs, flux, filt_files)
+        self.photometry_samples.append(phot)
 
         return phot
 
@@ -68,6 +70,7 @@ class galaxyfitter(object):
     # Fit an array of photometry, output the best galaxy parameters
     def fit(self, photdata, prior, filt_files):
 
+
         filt = filter_set(filt_files)
         effwavs = filt.eff_wavs
 
@@ -79,9 +82,13 @@ class galaxyfitter(object):
         dat[0] = 2.99792458E-05 * ((1e-6 * dat[0]) / ((effwavs)**2))
         dat[1] = 2.99792458E-05 * ((1e-6 * dat[1]) / ((effwavs)**2))
 
+
         # Run the nested sampler over the likelihood function
+        self.spectrum_samples=[]
+        self.photometry_samples=[]
         sampler = Sampler(prior, lambda param_dict: self.log_likelihood(dat, self.gen_phot, filt_files, param_dict), n_live=500)
-        success = sampler.run(verbose=True, discard_exploration=True, timeout=1800)
+        success = sampler.run(verbose=True, discard_exploration=False, timeout=1800)
+
 
         # Plot results
         points, log_w, log_l = sampler.posterior()
@@ -94,6 +101,11 @@ class galaxyfitter(object):
             samples = np.random.choice(points_i, size = points_i.shape[0], p=np.exp(log_w))
             dic[key + "_50"] = np.percentile(samples, 50)
 
+
+        choicer = np.random.choice(np.arange(len(self.spectrum_samples)), size = len(self.spectrum_samples), p=np.exp(log_w))
+        self.spectrum_samples = np.take(np.array((self.spectrum_samples)), choicer, axis=0)
+        self.photometry_samples = np.take(np.array((self.photometry_samples)), choicer, axis=0)
+
         exp = {}
         exp["age"] = dic['age_50']
         exp["tau"] = dic['tau_50']
@@ -106,10 +118,12 @@ class galaxyfitter(object):
         fit_instructions["exponential"] = exp   
         fit_instructions["dust"] = dust
 
+
         print('Finished in ' + str(time.time() - t0)+ ' seconds.')
         print(' ')
 
-        return fit_instructions
+
+        return fit_instructions, self.wavs, self.spectrum_samples, self.photometry_samples
 
 
 

{chanter-0.0.9 → chanter-0.1.0}/chanter/model_galaxy.py

@@ -56,7 +56,10 @@ igm_grid = resample_igm_grid(rest_target_wavs)
 
 
 '''
-# !!AGE RESAMPLING NOT WORKING!!
+######
+# Age resampling currently produces strange results. This is hence deactivated for now, which only serves to slow the code down and does
+# not affect accuracy of results.
+######
 
 def resample_ssp_ages(grid_raw_ages, target_ages):
 

chanter-0.1.0/chanter.egg-info/PKG-INFO

@@ -0,0 +1,25 @@
+Metadata-Version: 2.4
+Name: chanter
+Version: 0.1.0
+Summary: Simple quiescent galaxy spectral modelling and fitting
+Author: Struan Stevenson
+Author-email: struan.stevenson@ed.ac.uk
+Description-Content-Type: text/markdown
+Requires-Dist: numpy
+Requires-Dist: astropy
+Requires-Dist: matplotlib
+Requires-Dist: spectres
+Requires-Dist: nautilus-sampler
+Requires-Dist: json
+Dynamic: author
+Dynamic: author-email
+Dynamic: description
+Dynamic: description-content-type
+Dynamic: requires-dist
+Dynamic: summary
+
+# CHANTER
+
+#### CHANTER is a simple quiecent galaxy spectral energy distribution (SED) modelling/fitting code, based on the larger BAGPIPES code, of astronomical literary fame. CHANTER provides the barebones of Stellar Population Synthesis (SPS), as well as redshift and dust/IGM attenuation effects, and is written in a simple and digestible format, ideal for learning the theory behind larger spectral fitting software codes. CHANTER does not yet include nebular attenuation, and therefore only models the SED of quiescent galaxies.
+
+Currently, I cannot include the SSP/IGM grids in the package due to PyPi size limitations, but this is being appealed.

{chanter-0.0.9 → chanter-0.1.0}/chanter.egg-info/requires.txt

@@ -3,3 +3,4 @@ astropy
 matplotlib
 spectres
 nautilus-sampler
+json

{chanter-0.0.9 → chanter-0.1.0}/setup.py

@@ -8,9 +8,9 @@ long_description = (this_directory / "README.md").read_text()
 setup(
     name='chanter',
 
-    version='0.0.9',
+    version='0.1.0',
 
-    description='Simple galaxy spectral modelling and fitting',
+    description='Simple quiescent galaxy spectral modelling and fitting',
 
     long_description=long_description,
 
@@ -24,6 +24,6 @@ setup(
 
     package_data = {'': ['*.txt', '*.fits'],},  
 
-    install_requires=["numpy", "astropy", "matplotlib", "spectres", "nautilus-sampler"],
+    install_requires=["numpy", "astropy", "matplotlib", "spectres", "nautilus-sampler", "json"],
 
 )

chanter-0.0.9/PKG-INFO

@@ -1,24 +0,0 @@
-Metadata-Version: 2.4
-Name: chanter
-Version: 0.0.9
-Summary: Simple galaxy spectral modelling and fitting
-Author: Struan Stevenson
-Author-email: struan.stevenson@ed.ac.uk
-Description-Content-Type: text/markdown
-Requires-Dist: numpy
-Requires-Dist: astropy
-Requires-Dist: matplotlib
-Requires-Dist: spectres
-Requires-Dist: nautilus-sampler
-Dynamic: author
-Dynamic: author-email
-Dynamic: description
-Dynamic: description-content-type
-Dynamic: requires-dist
-Dynamic: summary
-
-# CHANTER
-
-#### CHANTER is a simple galaxy spectral energy distribution (SED) modelling/fitting code, based on the larger BAGPIPES code, of astronomical literary fame. CHANTER provides the barebones of Stellar Population Synthesis (SPS), as well as redshift and dust/IGM attenuation effects, and is written in a simple and digestible format, ideal for learning the theory behind larger spectral fitting software codes.
-
-Currently, I cannot include the SSP/IGM grids in the package due to PyPi size limitations, but this is being appealed.

chanter-0.0.9/README.md

@@ -1,5 +0,0 @@
-# CHANTER
-
-#### CHANTER is a simple galaxy spectral energy distribution (SED) modelling/fitting code, based on the larger BAGPIPES code, of astronomical literary fame. CHANTER provides the barebones of Stellar Population Synthesis (SPS), as well as redshift and dust/IGM attenuation effects, and is written in a simple and digestible format, ideal for learning the theory behind larger spectral fitting software codes.
-
-Currently, I cannot include the SSP/IGM grids in the package due to PyPi size limitations, but this is being appealed.

chanter-0.0.9/chanter.egg-info/PKG-INFO

@@ -1,24 +0,0 @@
-Metadata-Version: 2.4
-Name: chanter
-Version: 0.0.9
-Summary: Simple galaxy spectral modelling and fitting
-Author: Struan Stevenson
-Author-email: struan.stevenson@ed.ac.uk
-Description-Content-Type: text/markdown
-Requires-Dist: numpy
-Requires-Dist: astropy
-Requires-Dist: matplotlib
-Requires-Dist: spectres
-Requires-Dist: nautilus-sampler
-Dynamic: author
-Dynamic: author-email
-Dynamic: description
-Dynamic: description-content-type
-Dynamic: requires-dist
-Dynamic: summary
-
-# CHANTER
-
-#### CHANTER is a simple galaxy spectral energy distribution (SED) modelling/fitting code, based on the larger BAGPIPES code, of astronomical literary fame. CHANTER provides the barebones of Stellar Population Synthesis (SPS), as well as redshift and dust/IGM attenuation effects, and is written in a simple and digestible format, ideal for learning the theory behind larger spectral fitting software codes.
-
-Currently, I cannot include the SSP/IGM grids in the package due to PyPi size limitations, but this is being appealed.