| Title: | Evolutionary Tracks and Isochrones from Pisa Stellar Evolution Database |
|---|---|
| Description: | Manages and display stellar tracks and isochrones from Pisa low-mass database. Includes tools for isochrones construction and tracks interpolation. |
| Authors: | Matteo Dell'Omodarme [aut, cre], Giada Valle [aut] |
| Maintainer: | Matteo Dell'Omodarme <[email protected]> |
| License: | GPL (>= 2) |
| Version: | 0.3-4 |
| Built: | 2025-03-01 05:51:42 UTC |
| Source: | https://github.com/cran/stellaR |
The package is designed to manage and display stellar tracks and isochrones from Pisa low-mass database. Data are archived at CDS and are available for download (http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/540/A26).
The package includes tools to gather data from an active Internet connection, print and plot them with appropriate format.
Tools for make new isocrones at desired ages, and to interpolate the database for composition not already available are also provided.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
The data sets gives the compositions of the computed models.
compositionscompositions
A list containing 4 variables:
the initial metallicities of the models.
a data frame with the initial helium abundances of the models.
the mixing-lengths of the models.
the alpha-enhanchments of the models. It is coded
0 for [alpha/Fe] = 0.0 and 1 for
[alpha/Fe] = 0.3.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
Import stellar track data for a grid of HB models starting from ZAHB and ending at the onset of thermal pulses.
getHbgrid(z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")getHbgrid(z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")
z |
the initial metallicity of the model to get. Allowed values are 0.0001, 0.0002, 0.0003, 0.0004, 0.0005, 0.0006, 0.0007, 0.0008, 0.0009, 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009, 0.01. |
y |
the initial helium abundance of the model to get. Allowed values are 0.25, 0.27, 0.33, 0.38, 0.42. |
ml |
the mixing-length of the model to get. Allowed values are 1.7, 1.8, 1.9. |
afe |
the alpha-enhanchment of the model to get. It can be one of
|
baseURL |
the URL of the base directory of the database from where get the models. |
As a sanity check, the function performs a consistency test on the parameters
z, y, ml, afe by a call to
testComposition.
By default the function requires an Internet access for data
gathering. Tracks and isochrones are downloaded from the on-line
repository CDS
(http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/540/A26).
It is also possible to gather data from a local download of the
catalog, specifying as baseURL the root directory of the
download. In this case the structure of the catalog and the names of
files and directory must be preserved.
getHbgrid returns an object of class hbset, i.e.
a list with one element of class hb for each HB model.
If baseURL points to an invalid path or CDS ftp site is
unavailable
the function returns NA and issues a warning.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
getTrk, getHb, getIso, getTrkSet.
### slow! ## Not run: hbgrid <- getHbgrid(0.002, 0.25, 1.7, 0) ### get data from local directory /data hbgrid <- getHbgrid(0.002, 0.25, 1.7, 0, baseURL="/data/") ## End(Not run)### slow! ## Not run: hbgrid <- getHbgrid(0.002, 0.25, 1.7, 0) ### get data from local directory /data hbgrid <- getHbgrid(0.002, 0.25, 1.7, 0, baseURL="/data/") ## End(Not run)
Import a stellar isochrone.
getIso(age, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")getIso(age, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")
age |
the age (in Gyr) of the isochrone to get. Allowed values are 8.0, 8.5, 9.0, 9.5, 10.0, 10.5, 11.0, 11.5, 12.0, 12.5, 13.0, 13.5, 14.0, 14.5, 15.0. |
z |
the initial metallicity of the model to get. Allowed values are 0.0001, 0.0002, 0.0003, 0.0004, 0.0005, 0.0006, 0.0007, 0.0008, 0.0009, 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009, 0.01. |
y |
the initial helium abundance of the model to get. Allowed values are 0.25, 0.27, 0.33, 0.38, 0.42. |
ml |
the mixing-length of the model to get. Allowed values are 1.7, 1.8, 1.9. |
afe |
the alpha-enhanchment of the model to get. It can be one of
|
baseURL |
the URL of the base directory of the database from where get the models. |
As a sanity check, the function performs a consistency test on the parameters
z, y, ml, afe by a call to
testComposition.
By default the function requires an Internet access for data
gathering. Tracks and isochrones are downloaded from the on-line
repository CDS
(http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/540/A26).
It is also possible to gather data from a local download of the
catalog, specifying as baseURL the root directory of the
download. In this case the structure of the catalog and the names of
files and directory must be preserved.
getIso returns an object of class iso.
More specifically this is a list with elements:
age |
the age of the isochrone (in Gyr) |
z |
the metallicity of the model |
y |
the initial helium content of the model |
alpha.enh |
the mixing-length value of the model |
data |
the dataset of stellar isocrone, containing the following five variables:
|
print, and plot methods are available for
the class.
If baseURL points to an invalid path or CDS ftp site is
unavailable
the function returns NA and issues a warning.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
getZahb, getHbgrid, getTrk, getTrkSet.
## Not run: iso <- getIso(12.0, 0.002, 0.25, 1.7, 0) ### get data from local directory /data iso <- getIso(12.0, 0.002, 0.25, 1.7, 0, baseURL="/data/") ## End(Not run)## Not run: iso <- getIso(12.0, 0.002, 0.25, 1.7, 0) ### get data from local directory /data iso <- getIso(12.0, 0.002, 0.25, 1.7, 0, baseURL="/data/") ## End(Not run)
Import a stellar evolutionary track from pre-main sequence to He flash and from ZAHB to thermal pulses.
getTrk(m, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/") getHb(m, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")getTrk(m, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/") getHb(m, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")
m |
the mass, in unit of solar mass, of the model to get. Allowed values are 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00, 1.05, 1.10. |
z |
the initial metallicity of the model to get. Allowed values are 0.0001, 0.0002, 0.0003, 0.0004, 0.0005, 0.0006, 0.0007, 0.0008, 0.0009, 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009, 0.01. |
y |
the initial helium abundance of the model to get. Allowed values are 0.25, 0.27, 0.33, 0.38, 0.42. |
ml |
the mixing-length of the model to get. Allowed values are 1.7, 1.8, 1.9. |
afe |
the alpha-enhanchment of the model to get. It can be one of
|
baseURL |
the URL of the base directory of the database from where
get the models. See |
As a sanity check, the function performs a consistency test on the parameters
z, y, ml, afe by a call to
testComposition.
By default the function requires an Internet access for data
gathering. Tracks and isochrones are downloaded from the on-line
repository CDS
(http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/540/A26).
It is also possible to gather data from a local download of the
catalog, specifying as baseURL the root directory of the
download. In this case the structure of the catalog and the names of
files and directory must be preserved.
getTrk returns an object of class trk; getHb
returns an object of class hb.
More specifically these objects are lists with common elements:
mass |
the mass of the model (in unit of solar mass) |
z |
the metallicity of the model |
y |
the initial helium content of the model |
ml |
the mixing-length value of the model |
alpha.enh |
the alpha-enhanchment of the model. It can be one of
|
data |
the dataset of stellar track, containing the following variables:
Only |
hb class contains also the element:
massRGB |
the mass of the RGB progenitor |
print, and plot methods are available for
the two classes.
If baseURL points to an invalid path or CDS ftp site is
unavailable
the function returns NA and issues a warning.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
getZahb, getHbgrid, getIso, getTrkSet.
## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) ### get data from local directory /data trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0, baseURL="/data/") ### multi-panel plot of the various quantities with time track <- getTrk(0.80, 0.001, 0.25, 1.90, 0) if(!is.na(track)[1]) { trkdata <- within(stack(track$data), time <- rep(track$data$time, length.out=length(values)) ) require(lattice) xyplot( values ~ time | ind, data=trkdata, type="l", scales=list(y=list(relation="free"))) } ## End(Not run)## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) ### get data from local directory /data trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0, baseURL="/data/") ### multi-panel plot of the various quantities with time track <- getTrk(0.80, 0.001, 0.25, 1.90, 0) if(!is.na(track)[1]) { trkdata <- within(stack(track$data), time <- rep(track$data$time, length.out=length(values)) ) require(lattice) xyplot( values ~ time | ind, data=trkdata, type="l", scales=list(y=list(relation="free"))) } ## End(Not run)
Import a set of stellar evolutionary tracks from pre-main sequence to He flash, or a set of isochrones.
getTrkSet(m, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/") getIsoSet(age, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")getTrkSet(m, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/") getIsoSet(age, z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")
m |
a vector of masses, in unit of solar mass, of the models to get. Allowed values are 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00, 1.05, 1.10. |
age |
a vector of ages, in Gyr, of the models to get. Allowed values are 8.0, 8.5, 9.0, 9.5, 10.0, 10.5, 11.0, 11.5, 12.0, 12.5, 13.0, 13.5, 14.0, 14.5, 15.0. |
z |
a vector of the initial metallicity of the models to get. Allowed values are 0.0001, 0.0002, 0.0003, 0.0004, 0.0005, 0.0006, 0.0007, 0.0008, 0.0009, 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009, 0.01. |
y |
a vector of the initial helium abundance of the models to get. Allowed values are 0.25, 0.27, 0.33, 0.38, 0.42. |
ml |
a vector of the mixing-length of the models to get. Allowed values are 1.7, 1.8, 1.9. |
afe |
the alpha-enhanchment of the models to get. It can be one of
|
baseURL |
the URL of the base directory of the database from where get the models. |
As a sanity check, the function performs a consistency test on the parameters
z, y, ml, afe by a call to
testComposition.
By default the function requires an Internet access for data
gathering. Tracks and isochrones are downloaded from the on-line
repository CDS
(http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/540/A26).
It is also possible to gather data from a local download of the
catalog, specifying as baseURL the root directory of the
download. In this case the structure of the catalog and the names of
files and directory must be preserved.
getTrkSet returns an object of class trkset, i.e. a
list of objects of class trk.
getIsoSet returns an object of class isoset, i.e. a
list of objects of class iso.
print, and plot methods are available for
the two classes.
If baseURL points to an invalid path or CDS ftp site is
unavailable
the function returns NA and issues a warning.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
getZahb, getHbgrid,
getHb, getIso.
## Not run: ### get two masses trkset <- getTrkSet(c(0.9, 1.0), 0.002, 0.25, 1.7, 0) ### get two masses at two metallicity, for a total of 4 objects trkset <- getTrkSet(c(0.9, 1.0), c(0.002, 0.01), 0.25, 1.7, 0) ### get data from local directory /data set <- getTrkSet(0.9, 0.002, c(0.25, 0.33), 1.7, 0, baseURL="/data/") ## End(Not run)## Not run: ### get two masses trkset <- getTrkSet(c(0.9, 1.0), 0.002, 0.25, 1.7, 0) ### get two masses at two metallicity, for a total of 4 objects trkset <- getTrkSet(c(0.9, 1.0), c(0.002, 0.01), 0.25, 1.7, 0) ### get data from local directory /data set <- getTrkSet(0.9, 0.002, c(0.25, 0.33), 1.7, 0, baseURL="/data/") ## End(Not run)
Import a stellar ZAHB model.
getZahb(z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")getZahb(z, y, ml, afe, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")
z |
the initial metallicity of the model to get. Allowed values are 0.0001, 0.0002, 0.0003, 0.0004, 0.0005, 0.0006, 0.0007, 0.0008, 0.0009, 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009, 0.01. |
y |
the initial helium abundance of the model to get. Allowed values are 0.25, 0.27, 0.33, 0.38, 0.42. |
ml |
the mixing-length of the model to get. Allowed values are 1.7, 1.8, 1.9. |
afe |
the alpha-enhanchment of the model to get. It can be one of
|
baseURL |
the URL of the base directory of the database from where get the models. |
As a sanity check, the function performs a consistency test on the parameters
z, y, ml, afe by a call to
testComposition.
By default the function requires an Internet access for data
gathering. Tracks and isochrones are downloaded from the on-line
repository CDS
(http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/540/A26)
It is also possible to gather data from a local download of the
catalog, specifying as baseURL the root directory of the
download. In this case the structure of the catalog and the names of
files and directory must be preserved.
getZahb returns an object of class zahb
More specifically it is a list with elements:
z |
the metallicity of the model |
y |
the initial helium content of the model |
ml |
the mixing-length value of the model |
alpha.enh |
the alpha-enhanchment of the model. It can be one of
|
data |
the dataset of stellar track, containing the following variables:
|
print, and plot methods are available for
the two classes.
If baseURL points to an invalid path or CDS ftp site is
unavailable
the function returns NA and issues a warning.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
getTrk, getHbgrid, getIso, getTrkSet.
## Not run: zahb <- getZahb(0.002, 0.25, 1.7, 0) ### get data from local directory /data zahb <- getZahb(0.002, 0.25, 1.7, 0, baseURL="/data/") ## End(Not run)## Not run: zahb <- getZahb(0.002, 0.25, 1.7, 0) ### get data from local directory /data zahb <- getZahb(0.002, 0.25, 1.7, 0, baseURL="/data/") ## End(Not run)
Construct a set of interpolated tracks from the objects in the database.
interpTrk(z, y, ml, afe, vmass=seq(0.30,1.10, by=0.05), baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")interpTrk(z, y, ml, afe, vmass=seq(0.30,1.10, by=0.05), baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")
z |
the initial metallicity of the models to generate. Allowed values
are between |
y |
the initial helium abundance of the models to generate. Allowed
values are between |
ml |
the mixing-length of the model to generate. Allowed values are
between |
afe |
the alpha-enhanchment of the model to get. It can be one of
|
vmass |
vector of masses of the stellar model to generate. Default to the whole set in the database. Allowed values are 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, 0.80, 0.85, 0.90, 0.95, 1.00, 1.05, 1.10. |
baseURL |
the URL of the base directory of the database from where get the models. |
A 3D interpolation is performed on the database to construct the required set of tracks. Interpolation on helium abundance and on mixing-length occur linearly on the variables, while the metallicity is log-transformed before the interpolation procedure.
The function returns an object of class trkset containing the
interpolated tracks.
### slow! ## Not run: iptrk <- interpTrk(0.002, 0.25, 1.74, 0) ### get data from local directory /data iptrk <- interpTrk(0.002, 0.25, 1.74, 0, baseURL="/data/") ## End(Not run)### slow! ## Not run: iptrk <- interpTrk(0.002, 0.25, 1.74, 0) ### get data from local directory /data iptrk <- interpTrk(0.002, 0.25, 1.74, 0, baseURL="/data/") ## End(Not run)
Extract some relevant evolutionary points from stellar tracks or stellar isochrones.
## S3 method for class 'trk' keypoints(x, ...) ## S3 method for class 'trkset' keypoints(x, ...) ## S3 method for class 'iso' keypoints(x, ...) ## S3 method for class 'isoset' keypoints(x, ...)## S3 method for class 'trk' keypoints(x, ...) ## S3 method for class 'trkset' keypoints(x, ...) ## S3 method for class 'iso' keypoints(x, ...) ## S3 method for class 'isoset' keypoints(x, ...)
x |
an object of classes |
... |
further arguments, not implemented. |
keypoints.trk extracts from the object x$data the rows
relative to the following evolutionary stages:
ZAMS: Zero-Age Main-Sequence, defined as the point for which the central H abundance drops below 99% of its initial value.
TO: Turn-Off, defined as the point for which the effective temperature reaches its maximum value. If multiple lines satisfy the constraint, the values of all the rows are averaged.
BTO: Brighter Turn-Off, defined as the point for which the effective temperature drops below the one of the TO minus 100 K. The points could not exist for low masses. See Chaboyer et al. (1996) for an analysis of the advantages of this evolutionary point with respect to TO.
exHc: Central H exhaustion, defined as the point for which the central H abundance is zero. For low masses the point could coincide with TO. This is the last point of the tracks with mass lower or equal to 0.50 Msun.
Heflash: Helium flash, the last point of the track for masses higher than 0.50 Msun.
keypoints.trkset performs the selection described above for all the
set of objects.
The returned object contains the variable id, which labels the
evolutionary phases with the following coding: 1 = ZAMS, 2 = TO, 3 =
BTO, 4 = exHc, 5 = Heflash.
keypoints.iso extracts from the object x$data the rows
relative to the following evolutionary stages:
TO: Turn-Off, defined as the point for which the effective temperature reaches its maximum value. If multiple lines satisfy the constraint, the values of all the rows are averaged.
BTO: Brighter Turn-Off, defined as the point for which the effective temperature drops below the one of the TO minus 100 K.
keypoints.isoset performs the selection described above for all the
set of objects.
The returned object contains the variable id, which labels the
evolutionary phases with the following coding: 1 = TO, 2 = BTO.
A data frame containing the observations relative to the selected evolutionary phases. The columns relative to mass (or age for isochrones), metallicity, initial He value, mixing-length, alpha-enhancement, and phase identifier are inserted in the data frame.
B. Chaboyer, P. Demarque, P.J. Kernan, L.M. Krauss, and A. Sarajedini (1996). An accurate relative age estimator for globular clusters. MNRAS, 283, 683-689.
## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(trk)[1]) keypoints(trk) is <- getIso(11, 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(is)[1]) keypoints(is) ## End(Not run)## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(trk)[1]) keypoints(trk) is <- getIso(11, 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(is)[1]) keypoints(is) ## End(Not run)
The function computes by interpolation stellar isochrones of given ages from evolutionary tracks.
makeIso(age, z=NULL, y=NULL, ml=NULL, afe=NULL, log=FALSE, linear=TRUE, tr=NULL, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")makeIso(age, z=NULL, y=NULL, ml=NULL, afe=NULL, log=FALSE, linear=TRUE, tr=NULL, baseURL="ftp://cdsarc.u-strasbg.fr/pub/cats/J/A+A/540/A26/")
age |
the ages (in Gyr) of the isochrones to construct. A vector of ages is allowed. |
z |
the initial metallicity of the isochrones to build. |
y |
the initial helium abundance of the isochrones to build. |
ml |
the mixing-length of the isochrones to build. |
afe |
the alpha-enhanchment of the isochrone to build. It can be one of
|
log |
a logical value indicating whether the interpolation of the tracks should be performed on the logarithm of the evolutionary time. |
linear |
a logical value indicating whether linear interpolation
should be performed. Only the option |
tr |
a set of track to be used for isocrones construction. See
|
baseURL |
the URL of the base directory of the database from where get the models. |
Isochrones are obtained by mean of a standard interpolation procedure on
the set of tracks. Let be the set of tracks, parametrized
by the value of the mass .
Let be the evolutionary
time for the ith point on the track of mass .
Let be the point on the track of lower mass of for
which is greater of the time required for the isochrone.
For each point on , an interpolation of mass,
logarithm of the effective temperature and logarithm of the luminosity
is performed among tracks. These points define the required isochrone.
If a set of tracks is supplied by mean of the argument tr, the
function interpolates among these tracks. The values of
z, y, ml, afe are recovered from the
supplied objects and a test of consistency is performed to assure that
the tracks are homogeneous in these parameters.
The function returns an object of class isoset.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
### slow! ## Not run: isoset <- makeIso(11.2, 0.002, 0.25, 1.7, 0) ## End(Not run) ### get data from local directory /data ## Not run: isoset <- makeIso(c(11.2, 12.4), 0.002, 0.25, 1.7, 0, baseURL="/data/")### slow! ## Not run: isoset <- makeIso(11.2, 0.002, 0.25, 1.7, 0) ## End(Not run) ### get data from local directory /data ## Not run: isoset <- makeIso(c(11.2, 12.4), 0.002, 0.25, 1.7, 0, baseURL="/data/")
These data sets give the stellar mass in HB for the progenitors in RGB.
masshb masshbgridmasshb masshbgrid
The two data frames contain 6 variables:
the RGB mass of the model.
the initial metallicity of the model.
the initial helium abundance of the model.
the mixing-length of the model.
the alpha-enhanchment of the model. It is coded
mix=AS09a0 for [alpha/Fe] = 0.0 and mix=AS09a3 for
[alpha/Fe] = 0.3.
the HB mass of the model. The dataset
masshbgrid contains multiple values for masshb for
each single RGB progenitor mass.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
Plot one evolutionary track or one isochrone extracted from objects of
classes trk,
hb, zahb, iso.
## S3 method for class 'trk' plot(x, add=FALSE, ...) ## S3 method for class 'hb' plot(x, add=FALSE, ...) ## S3 method for class 'zahb' plot(x, add=FALSE, ...) ## S3 method for class 'iso' plot(x, add=FALSE, ...)## S3 method for class 'trk' plot(x, add=FALSE, ...) ## S3 method for class 'hb' plot(x, add=FALSE, ...) ## S3 method for class 'zahb' plot(x, add=FALSE, ...) ## S3 method for class 'iso' plot(x, add=FALSE, ...)
x |
an object of classes |
add |
a logical value indicating whether the plot should be added to an existing device. |
... |
further arguments passed to |
The contents of the vectors x$data$logTe and x$data$logL are
plotted by a call to plotAstro. Sensible values for the
axes major ticks and labels are chosen.
## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(trk)[1]) plot(trk) ## End(Not run)## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(trk)[1]) plot(trk) ## End(Not run)
Plot a set of evolutionary tracks from objects of
classes trkset, and
hbset, or a set of isochrones from object of class isoset.
## S3 method for class 'trkset' plot(x, add=FALSE, col=1, lty=1, xlim=NULL, ylim=NULL, ...) ## S3 method for class 'hbset' plot(x, add=FALSE, col=1, lty=1, xlim=NULL, ylim=NULL, ...) ## S3 method for class 'isoset' plot(x, add=FALSE, col=1, lty=1, xlim=NULL, ylim=NULL, ...)## S3 method for class 'trkset' plot(x, add=FALSE, col=1, lty=1, xlim=NULL, ylim=NULL, ...) ## S3 method for class 'hbset' plot(x, add=FALSE, col=1, lty=1, xlim=NULL, ylim=NULL, ...) ## S3 method for class 'isoset' plot(x, add=FALSE, col=1, lty=1, xlim=NULL, ylim=NULL, ...)
x |
an object of class |
, or isoset.
add |
a logical value indicating whether the plot should be added to an existing device. |
col |
the colors for lines. Multiple colors can be specified so that each track can be given its own color. If there are fewer colors than tracks they are recycled in the standard fashion. |
lty |
the line types. Multiple types can be specified so that each track can be given its own type. If there are fewer type than tracks they are recycled in the standard fashion. |
xlim |
the x limits (x1, x2) of the plot.
The default value, |
ylim |
the y limits of the plot. |
... |
further arguments passed to |
For each object X contained in the input of class trkset,
hbset, or
isoset the contents
of the vectors X$data$logTe and X$data$logL are
plotted by a call to plotAstro. Sensible values for the
axes ranges, major ticks and labels are chosen.
## Not run: trkset <- getTrkSet(c(0.7, 0.9, 1.1), 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(trkset)[1]) plot(trkset) ## End(Not run)## Not run: trkset <- getTrkSet(c(0.7, 0.9, 1.1), 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(trkset)[1]) plot(trkset) ## End(Not run)
Draw a scatter plot in the active graphics window.
plotAstro(x, y, type="l", xlab="", ylab="", xi=(max(x)-min(x))/5, yi=(max(y)-min(y))/5, xmt=3, ymt=3, revX=FALSE, revY=FALSE, xlim=NULL, ylim=NULL, cex=1.0, cex.axis=1.3, cex.lab=1.5, add=FALSE, ...)plotAstro(x, y, type="l", xlab="", ylab="", xi=(max(x)-min(x))/5, yi=(max(y)-min(y))/5, xmt=3, ymt=3, revX=FALSE, revY=FALSE, xlim=NULL, ylim=NULL, cex=1.0, cex.axis=1.3, cex.lab=1.5, add=FALSE, ...)
x, y
|
the |
type |
1-character string giving the type of plot desired. The
following values are possible: |
xlab |
a label for the x axis, default to empty string. |
ylab |
a label for the y axis, default to empty string. |
xi |
major ticks spacing for x axis. |
yi |
major ticks spacing for y axis. |
xmt |
number of minor ticks between two major ticks for x axis. |
ymt |
number of minor ticks between two major ticks for y axis. |
revX |
a logical value indicating whether the x axis should be reversed. |
revY |
a logical value indicating whether the y axis should be reversed. |
xlim |
the x limits (x1, x2) of the plot.
The default value, |
ylim |
the y limits of the plot. |
cex |
the amount by which plotting text and symbols should be magnified relative to the default. |
cex.axis |
the magnification to be used for axis annotation
relative to the current setting of |
cex.lab |
the magnification to be used for x and y labels
relative to the current setting of |
add |
a logical value indicating whether the plot should be added to an existing device. |
... |
further arguments passed to low level plot functions. |
## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(trk)[1]) plot(trk) ## End(Not run)## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) ### check return value from CDS if(!is.na(trk)[1]) plot(trk) ## End(Not run)
Format and print the contents of the objects of classes trk,
hb, zahb, iso.
## S3 method for class 'trk' print(x, ...) ## S3 method for class 'hb' print(x, ...) ## S3 method for class 'zahb' print(x, ...) ## S3 method for class 'iso' print(x, ...)## S3 method for class 'trk' print(x, ...) ## S3 method for class 'hb' print(x, ...) ## S3 method for class 'zahb' print(x, ...) ## S3 method for class 'iso' print(x, ...)
x |
An object of classes |
... |
further arguments passed to low level print functions. |
The contents of the variables mass (or age), z, y,
ml, alpha.enh (and massRGB if appliable)
from classes trk,
hb, zahb, or iso
are fomatted and printed. The function returns its argument invisibly (via invisible(x)).
## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) trk ## End(Not run)## Not run: trk <- getTrk(0.9, 0.002, 0.25, 1.7, 0) trk ## End(Not run)
Show the values of chemical composition (Z, Y, AFE) and mixing-length present in the database.
showComposition()showComposition()
Print the combinations of the inputs that exist in the database.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
showComposition()showComposition()
Check that the given chemical composition (Z, Y, [alpha/Fe]) and the mixing-length value exist in the stellar database.
testComposition(Z, Y, ML, AFE)testComposition(Z, Y, ML, AFE)
Z |
the initial metallicity of the model. |
Y |
the initial helium abundance of the model. |
ML |
the mixing-length of the model. |
AFE |
the alpha-enhanchment of the model. Allowed values are |
The function returns TRUE if the combination of the inputs
exists in the database.
M. Dell'Omodarme, G. Valle, S. Degl'Innocenti, and P.G. Prada Moroni (2012). The Pisa Stellar Evolution Data Base for low-mass stars. Astronomy and Astrophysics, 540, A26.
testComposition(0.002, 0.25, 1.7, 0)testComposition(0.002, 0.25, 1.7, 0)