Consistency Tests¶
Circle Test¶
The first test is to make sure the MCMC can get out the same result that FSPS was given. Seeing the flow of the values
| Input 1 (to FSPS) | Output 1 (from FSPS) | Input 2 (to spae) | Output 2 (from spea) |
|---|---|---|---|
| SFH parameters, redshift | SED | SED, redshift | SFH parameters |
So the output of spea should be the same as the input to FSPS.
Setting and generating the values in data/circlePhotometry.tsv¶
Lets generate the needed stellar population. This should be the same as the default in calcualteAge.pyL427.
import fsps
sp = fsps.StellarPopulation(zcontinuous=2, cloudy_dust=True,
add_neb_emission = True, sfh=5)
sdss_bands = ['sdss_u', 'sdss_g', 'sdss_r', 'sdss_i', 'sdss_z']
Here are the different input 1, star formation histories in a table
| name | logzsol | dust2 | tau | tStart | sfTrans | sfSlope | description |
|---|---|---|---|---|---|---|---|
| c1 | -0.5 | 0.1 | 0.5 | 1.5 | 9.0 | -1.0 | old |
| c2 | -0.5 | 0.1 | 0.5 | 1.5 | 9.0 | 15.0 | sharp burst & young |
| c3 | -0.5 | 0.1 | 7.0 | 3.0 | 10.0 | 15.0 | flat burst & young |
| c4 | -0.5 | 0.1 | 7.0 | 3.0 | 13.0 | 0.0 | flat burst |
| c5 | -1.5 | 0.1 | 0.5 | 1.5 | 9.0 | -1.0 | metal poor, like c1 |
| c6 | -0.5 | 0.8 | 7.0 | 3.0 | 10.0 | 15.0 | dusty, like c2 |
| c7 | -0.5 | 0.1 | 0.5 | 1.5 | 6.0 | 15.0 | mostly late linear, like c2 |
| c8 | -0.5 | 0.1 | 0.1 | 8.0 | 12.0 | 20.0 | very young |
and as copyable code
logzsol, dust2, tau, tStart, sfTrans, sfSlope = -0.5, 0.1, 0.5, 1.5, 9.0, -1.0
logzsol, dust2, tau, tStart, sfTrans, sfSlope = -0.5, 0.1, 0.5, 1.5, 9.0, 15.0
logzsol, dust2, tau, tStart, sfTrans, sfSlope = -0.5, 0.1, 7.0, 3.0, 10.0, 15.0
logzsol, dust2, tau, tStart, sfTrans, sfSlope = -0.5, 0.1, 7.0, 3.0, 13.0, 0.0
logzsol, dust2, tau, tStart, sfTrans, sfSlope = -1.5, 0.1, 0.5, 1.5, 9.0, -1.0
logzsol, dust2, tau, tStart, sfTrans, sfSlope = -0.5, 0.8, 7.0, 3.0, 10.0, 15.0
logzsol, dust2, tau, tStart, sfTrans, sfSlope = -0.5, 0.1, 0.5, 1.5, 6.0, 15.0
logzsol, dust2, tau, tStart, sfTrans, sfSlope = -0.5, 0.1, 0.1, 8.0, 12.0, 20.0
sp.params['logzsol'] = logzsol
dust1 = 2.0*dust2
sp.params['dust1'] = dust1
sp.params['dust2'] = dust2
sp.params['tau'] = tau
sp.params['sf_start'] = tStart
sp.params['sf_trunc'] = sfTrans
sp.params['sf_slope'] = sfSlope
sp.get_mags(tage=13.185, redshift=0.05, bands=sdss_bands)
And the results can be found in data/circlePhotometry.tsv available on GitHub with a boost of c = 25
comment on photometric uncertainties.
SED Check¶
Messier Objects¶
The paper coves this well.
Gupta Objects¶
The paper covers this well.