The SKIRT project
advanced radiative transfer for astrophysics
MarastonSEDFamily Class Reference

#include <MarastonSEDFamily.hpp>

Inheritance diagram for MarastonSEDFamily:

Public Types

enum class  IMF : int { Kroupa , Salpeter }

Public Member Functions

 MarastonSEDFamily (SimulationItem *parent, IMF imf)
double cdf (Array &lambdav, Array &pv, Array &Pv, const Range &wavelengthRange, const Array &parameters) const override
IMF imf () const
Range intrinsicWavelengthRange () const override
vector< SnapshotParameterparameterInfo () const override
double specificLuminosity (double wavelength, const Array &parameters) const override
Public Member Functions inherited from SimulationItem
template<class T>
T * find (bool setup=true) const
template<class T>
T * interface (int levels=-999999, bool setup=true) const
virtual string itemName () const
void setup ()
string typeAndName () const
Public Member Functions inherited from Item
 Item (const Item &)=delete
virtual ~Item ()
void addChild (Item *child)
const vector< Item * > & children () const
virtual void clearItemListProperty (const PropertyDef *property)
void destroyChild (Item *child)
virtual bool getBoolProperty (const PropertyDef *property) const
virtual vector< double > getDoubleListProperty (const PropertyDef *property) const
virtual double getDoubleProperty (const PropertyDef *property) const
virtual string getEnumProperty (const PropertyDef *property) const
virtual int getIntProperty (const PropertyDef *property) const
virtual vector< Item * > getItemListProperty (const PropertyDef *property) const
virtual ItemgetItemProperty (const PropertyDef *property) const
virtual string getStringProperty (const PropertyDef *property) const
int getUtilityProperty (string name) const
virtual void insertIntoItemListProperty (const PropertyDef *property, int index, Item *item)
Itemoperator= (const Item &)=delete
Itemparent () const
virtual void removeFromItemListProperty (const PropertyDef *property, int index)
virtual void setBoolProperty (const PropertyDef *property, bool value)
virtual void setDoubleListProperty (const PropertyDef *property, vector< double > value)
virtual void setDoubleProperty (const PropertyDef *property, double value)
virtual void setEnumProperty (const PropertyDef *property, string value)
virtual void setIntProperty (const PropertyDef *property, int value)
virtual void setItemProperty (const PropertyDef *property, Item *item)
virtual void setStringProperty (const PropertyDef *property, string value)
void setUtilityProperty (string name, int value)
virtual string type () const

Protected Member Functions

 MarastonSEDFamily ()
void setupSelfBefore () override
Protected Member Functions inherited from SEDFamily
 SEDFamily ()
Protected Member Functions inherited from SimulationItem
 SimulationItem ()
virtual bool offersInterface (const std::type_info &interfaceTypeInfo) const
virtual void setupSelfAfter ()
Protected Member Functions inherited from Item
 Item ()

Private Types

using BaseType
using ItemType

Private Attributes

IMF _imf
StoredTable< 3 > _table

Friends

class ItemRegistry

Detailed Description

An instance of the MarastonSEDFamily class represents the family of Maraston SEDs for single stellar populations (SSPs), parameterized on metallicity and age and scaled by the initial mass of the SSP (Maraston 1998, MNRAS, 300, 872–892). We use both the Kroupa and Salpeter models, in each case for the red horizontal branch morphology (the models for the blue horizontal branch morphology cover only a fraction of the parameter space). In other words, the SED family is available for two assumed initial mass functions (Kroupa and Salpeter).

For metallicities [Z/H] = -1.35, -0.33, 0, 0.35, models are given for ages from 10^3 yr to 15 Gyr. For metallicities [Z/H] = -2.25 and +0.67, only models from 1 to 15 Gyr are provided. These six metallicity values correspond to \(Z \approx 0.00011247, 0.00089337, 0.0093547, 0.02, 0.04477442, 0.09354703\).

The SEDs are tabulated over a wavelength range from 0.009 \(\mu\mathrm{m}\) to 160 \(\mu\mathrm{m}\) with the spectral resolution shown in the figure below.

The data were downloaded from http://www.icg.port.ac.uk/~maraston/Claudia's_Stellar_Population_Model.html and converted to SKIRT stored table format for inclusion as a resource file. The stored table is opened during setup, and it is subsequently interpolated to the desired parameters and wavelength grid when needed.

When imported from a text column file, the parameters for this SED family must appear in the following order in the specified default units (unless these units are overridden by column header info):

\[ M_\mathrm{init}\,(\mathrm{M}_\odot) \quad Z\,(\mathrm{dimensionless}) \quad t\,(\mathrm{yr}) \]

Member Enumeration Documentation

◆ IMF

enum class MarastonSEDFamily::IMF : int
strong

The enumeration type indicating the assumed initial mass function (IMF).

Kroupa : "Kroupa IMF".

Salpeter : "Salpeter IMF".

Constructor & Destructor Documentation

◆ MarastonSEDFamily() [1/2]

MarastonSEDFamily::MarastonSEDFamily ( )
inlineprotected

Default constructor for concrete Item subclass MarastonSEDFamily: "a Maraston SED family for single stellar populations".

◆ MarastonSEDFamily() [2/2]

MarastonSEDFamily::MarastonSEDFamily ( SimulationItem * parent,
IMF imf )
explicit

This constructor can be invoked programmatically by classes that use a hard-coded SED family (as opposed to selected through the ski file). Before the constructor returns, the newly created object is hooked up as a child to the specified parent in the simulation hierarchy (so it will automatically be deleted), and its setup() function has been called.

Member Function Documentation

◆ cdf()

double MarastonSEDFamily::cdf ( Array & lambdav,
Array & pv,
Array & Pv,
const Range & wavelengthRange,
const Array & parameters ) const
overridevirtual

This function constructs both the normalized probability density function (pdf) and the corresponding normalized cumulative distribution function (cdf) for the SED with the specified parameters over the specified wavelength range. The function returns the normalization factor. The number and type of parameters must match the information returned by the parameterInfo() function; if not the behavior is undefined.

Implements SEDFamily.

◆ imf()

IMF MarastonSEDFamily::imf ( ) const
inline

This function returns the value of the discoverable IMF enumeration property imf: "the assumed initial mass function".

The default value for this property is given by the conditional value expression "Kroupa".

◆ intrinsicWavelengthRange()

Range MarastonSEDFamily::intrinsicWavelengthRange ( ) const
overridevirtual

This function returns the intrinsic wavelength range of the SED family. It retrieves this range from the underlying stored table.

Implements SEDFamily.

◆ parameterInfo()

vector< SnapshotParameter > MarastonSEDFamily::parameterInfo ( ) const
overridevirtual

This function returns the number and type of parameters used by this particular SED family as a list of SnapshotParameter objects. Each of these objects specifies unit information and a human-readable descripton for the parameter.

Implements SEDFamily.

◆ setupSelfBefore()

void MarastonSEDFamily::setupSelfBefore ( )
overrideprotectedvirtual

This function opens the appropriate resource file (in SKIRT stored table format).

Reimplemented from SimulationItem.

◆ specificLuminosity()

double MarastonSEDFamily::specificLuminosity ( double wavelength,
const Array & parameters ) const
overridevirtual

This function returns the specific luminosity \(L_\lambda\) (i.e. radiative power per unit of wavelength) for the SED with the specified parameters at the specified wavelength, or zero if the wavelength is outside of the SED's intrinsic wavelength range. The number and type of parameters must match the information returned by the parameterInfo() function; if not the behavior is undefined.

Implements SEDFamily.


The documentation for this class was generated from the following file: