Frankel2018SFH#

class cogsworth.sfh.Frankel2018SFH(scale_length=None, scale_height=None, tsfr=<Quantity 6.8 Gyr>, alpha=0.3, Fm=-1, gradient=<Quantity -0.075 1 / kpc>, Rnow=<Quantity 8.7 kpc>, gamma=0.3, zsun=0.0142, galaxy_age=<Quantity 12. Gyr>, **kwargs)[source]#

Bases: StarFormationHistory

A star formation history for a component of the Milky Way, based on Frankel+2018.

Parameters are the same as StarFormationHistory but additionally with the following:

Parameters:
scale_heightQuantity [length], optional

Scale height of the disc, by default 0.3*u.kpc

tsfrQuantity [time], optional

Star formation timescale, by default 6.8*u.Gyr

alphafloat, optional

Disc inside-out growth parameter, by default 0.3

Fmint, optional

Metallicity at centre of disc at tm, by default -1

gradientQuantity [1/length], optional

Metallicity gradient, by default -0.075/u.kpc

RnowQuantity [length], optional

Radius at which present day metallicity is solar, by default 8.7*u.kpc

gammafloat, optional

Time dependence of chemical enrichment, by default 0.3

zsunfloat, optional

Solar metallicity, by default 0.0142

galaxy_ageQuantity [time], optional

Maximum lookback time, by default 12*u.Gyr

Methods Summary

draw_heights(size)

Draw heights from an exponential distribution with scale height given by the class attribute.

draw_phi(size)

draw_radii(size)

Inverse CDF sampling of galactocentric radii using Frankel+2018 Eq.

get_metallicity()

Methods Documentation

draw_heights(size)[source]#

Draw heights from an exponential distribution with scale height given by the class attribute.

Parameters:
sizeint

How many heights to draw

Returns:
zQuantity [length]

Random heights above the plane

draw_phi(size)[source]#
draw_radii(size)[source]#

Inverse CDF sampling of galactocentric radii using Frankel+2018 Eq. 5. The scale length is calculated using Eq. 6.

Parameters:
sizeint

How many radii to draw

Returns:
rhoQuantity [length]

Random Galactocentric radius

get_metallicity()[source]#