Reduction procedure
Since the first version of the library, the reduction procedure
described
in 2001A&A...369.1048P has been improved in several aspects but the
general philosophy remains the same. The basic steps in the reduction
are
to
(1) correct the different orders in the spectra for the blaze effect
and connect the orders together
(2) mask the telluric lines and spikes due to cosmic rays : in the
new version, the mask of the spikes due to cosmic rays has been
suppressed because it was found to alter real features. Because of it,
in the previous version of the library Lick indices measured on our
spectra presented a slope compared to measurements on Jones spectra,
the strongest features were found weaker than real (thanks to G.
Worthey who stressed our attention on this important problem).
(3) make the flux calibration
(4) build an interpolator which in turn is used to generate the grid
that feeds the population synthesis program (PEGASE.HR).
The strengthening of the spectra have been improved: artifacts near
Balmer's lines were corrected. The wavelength interval has been
extended
by adding 5 orders in the blue. It was 410-680 nm, it is now 400-680.
The correction for telluric lines has been improved. We corrected some
small bugs in the FITS keywords (in particular an unfortunate bad
rounding
of CRVAL1).
Flux calibration:
The main improvement in the flux calibration is due to the increase
of the number of external and internal comparisons. We made 2 observing
runs to have more intercomparisons. We have also improved the color
equations. Finally, we have improved the determination of the Galactic
extinction. It is determined from either Klare and Neckel, Hakkila et
al. or Chen et al. after a consistency test between these various 3D
model. Correction of the galactic extinction is primordial for
population
synthesis and quite critical because there exists a correlation in the
library between the spectral type and the extinction (both hot
stars and red giants are often in the disk and have high extinction).
Interpolator:
The main goal of project was to make a stellar library for synthesis
of stellar population. This require to provide a spectrum at any
desired
point of the parameter space, and actually the spectral systhesis tools
we are using, PEGASE.HR, uses a grid of spectra with a defined mesh.
The interpolator produces a spectrum at given Teff, log g and [Fe/H]
by evaluating polynomials in a way quite similar to the Lick fitting
functions
(Worthey 1994). The polynomials include with up to 21 terms and up to
the
power 4 in each of the parameters with a selection of cross terms. Each
wavelength element is represented with an independent polynomial but
the
functional form is the same for the whole spectrum. As well as for the
Lick fitting functions the parameter space is divided in a few
overlapping
regions where different polynomials are fitted. The overlap between
these
regions is used to give a smooth transition between the solutions
obtained
in each region. In the overlap the fluxes are a linear interpolation
between
the polynomials of each region. The terms choosen, and the validity
range
of the polynomials were chosen empirically by trial and errors.
The quality of the interpolator can be checked by comparing the
interpolated
fluxes with the observed fluxes for the spectra in the library. So, at
each iteration we searched these residuals for discrepancies and
systematic
effects and we modified the form of the polynomials and the validity
regions
in order to improve the model. The solution adopted in the present
version
reflect the intrinsic variation of the lines and the actual
distribution
of the stars in the space of parameters.
Because the library contains a large number of F-K stars, we have
also
weighted the spectra in order to re-equilibrate the distributions. The
polynomials were fitted by least-square with a kappa-sigma rejection of
outliers (un-corrected bad pixels in some spectra). In the new version
we have much more hot stars (>10000K) and cool stars (<4500K)
than in
previous version and the interpolator is therefore considerably
improved
with respect to the previous version.
Absolute and internal interpolators: The
interpolator
was first computed using the stars with known atmospheric parameters,
giving
each stars a weight function of the quality of these determination.
This
version is called the "absolute interpolator". Inverting the
interpolator
for a spectrum corresponding to a star with unknown parameters allows
to
determine these parameters. For the stars used to build the
interpolator,
the comparison between the parameters resulting for inversion and those
from our catalogue allow to detect errors on the input parameters of
some
stars (or errors on the identification of some stars at the telescope).
In some cases these errors could be corrected and the interpolator were
re-computed. This comparison also allowed to assess the weaknesses of
the
interpolator and drived its improvement (see above) and finally the
statistics
on this comparison reflect the errors on the parameters in our
catalogue
plus the errors due to the processing and modelling.
We then determined the parameters for all the stars using the
absolute
interpolator and used them to compute a new interpolator: the internal
interpolator. This internal interpolator may suffer from some biases
introduced
by the absolute interpolator, but it has smaller ramdom errors because
the errors on the input atmospheric parameters are decreased. It is
also
more stable because it is computed on more spectra, which is
particularly
sensitive in the less populated regions of the space of parameters. We
used the internal interpolator to compute the internal determination of
the atmospheric parameters listed in table 1.
Physical flux interpolator, correction for Galactic
extinction:
The interpolator was yet computed using the fluxes normalized to the
continuum,
because it is better suited to the determination of atmospheric
parameters.
The reason is that our flux calibration is not absolutely precise and
since
the information on the atmospheric parameters is essentially contained
in the lines, we wanted to avoid to bias the atmospheric parameters
with
errors in the flux calibration. But to make the population synthesis we
need to generate spectra calibrated in fluxes as emitted by the stars
(ie.
above the earth atmosphere and corrected for interstellar extinction).
The next step was therefore to correct galactic extinction. We used for
that the Schild (1977) extinction law (the choice of the extinction law
is not critical in this region of wavelengths) scaled with the E(B-V)
excess
listed in our catalogue. We build an interpolator for source flux
spectra
using the internal determination of the atmospheric parameters. Again
it
was possible to compare the results of the interpolation with all the
input
spectra. This allowed to detect spectra with discrepant SED, either due
to flux calibration, to erroneous value of the Galactic extinction. The
most discrepant case were closely examined, searching the literature
for
possible explanations or hints for errors on the parameters and
checking
the consistency of the estimates of the Galactic extinction and the
various
photometric measurements for the objects. Errors in the catalogues were
detected and corrected and in other cases we substituted the catalogued
value of the Galactic extinction with a value estimated from the
photometry
(and atmospheric parameters).
A final test on the quality of the source flux internal interpolator
was to compare the colour of the stars in different regions of the HR
diagram
with the colors measured on the spectra (Table xxx). The result of the
comparison is quite accpetable but we still see a systematic effect
that
we believe is due to the color equations used to convert between the
colors
measured on Elodie spectra, the Johnson colors and the Tycho colors. We
would like to solve this discrepancy in the near future.
Consistency and external comparisons:
Measurements
of Lick indices and determination of the atmospheric parameters are
described
in the article and the results are presented in.
Description of "calibration" files
These files, whose archive number starts with "99", are primarly
intended
to be used by the specialized reduction programs and by the pipeline
processing.
Their description may not be accurate enough for a standalone usage,
and
for such applications the users are invited to contact us for
additional
help.
LL_ELODIE/99001
Blaze function used as reference This is a "s2d" file produced by the
ELODIE
reduction program (TACOS). "s2d" means: spectra extracted by TACOS,
each
line is one order. The Blaze function is obtained from a Tungsten lamp
spectrum. Each order is smoothed using a degree 17 polynomial and is
normalized
to its maximum value.
LL_ELODIE/99002
Correction to the TACOS blaze (LL_ELODIE/99001) Dividing a spectrum by
a TACOS blaze produces a stair case result (consecutive orders do not
connect
because of the independent nomalizations). In order to connect the
orders
and to reduce the instrumental response to a "typical" stellar SED, we
computed this correction using 3 metal-deficient stars observed during
the same observing run. The original "s2d" spectra are divided by
LL_ELODIE/99001
and the continuum of each order is modelled with a degree 2 polynomial.
LL_ELODIE/99003
Pseudo-continuum mask It is computed after the whole archive has been
reduced.
In each continuum-normalized spectrum, a simple threshold is used to
diagnostic
if a given wavelength point belong to the pseudo-continuum. We choosed
0.97 as threshold. The corresonding wavelength point is given the value
1 or 0 if it belongs respectively to the continuum or not.
LL_ELODIE/99003
is the average of this mask over the whole archive. In some sense it
indicates
the probability for each wavelength point to be a continuum point. It
is
used as a weight for iterrating the computation of LL_ELODIE/99002 and
of the whole archive. (determining the continuum is required to fit and
correct variations of the blaze, see the text of the article).
LL_ELODIE/99004
Primary flux calibration. The first axis is the wavelength direction
The
first line is the instrumental response, The second the airmass
dependence
and the third the "haze" dependence (see text). LL_ELODIE/99004 is
obtained
after a comparison between the instrumental fluxes and templates from
the
datasets L1985BURN L1998GLUS and L1997ALEK
LL_ELODIE/99005
Smoothed flux calibration file Each line of LL_ELODIE/99004 is modeled
by a degree 17 polynomial in order to discard the wavy residuals due to
an approximate modeling of the spectral resolution of the templates
LL_ELODIE/99006
Secondary flux calibration. LL_ELODIE/99005 with correction after
pairwyse
comparisons
LL_ELODIE/99007
Corrected flux calibration file The spectra calibrated in physical flux
with LL_ELODIE/99005 are compared with spectra from L1984JACO L1996SERO
L96111KP1. The details of these comparison are given in
LL_ELODIE/99007.
The mean comparison is smoothed and used to compute a correction to the
instrumental response. LL_ELODIE/99006 serves to compute the final flux
calibration relations associated to each spectrum in the archive.
LL_ELODIE/99008
Comparison with other templates. The first axis is the wavelength
direction.
The 9 first lines contains respectively the comparisons with (1)
L96111KP1
and (2) L96111KP2 (Jone's library, respectively red and blue bands),
(3)
L1996SERO (OHP observations), (4) L1996SERO (CFH observations), (5)
L1984JACO
(Jacoby's library) (6) L1987KIEH (Kiehling 1987) (7) L1994DANK (Danks
\&
Dennefeld 1994) (8) L1983GUNN (Gunn & Stryker 1983 ) and (9)
L1984NSOA
(Solar flux atlas). The 10th line combines the first 9 lines. And the
11th
is the final correction obtained by smoothing the 10th.
LL_ELODIE/99010
Absolute TGM model The first axis is the wavelength direction. The
second
axis contains the 20 parameters of the TGM model fitted to the spectra
for which the atmosphric parameters are known. The functions of the
models
are: 1 Cst 2 T 3 Z 4 G 5 T**2 6 T**3 7 T**4 8 T*Z 9 T*G 10 G**2 11 Z**2
12 T**2*G 13 T**2*Z 14 G**3 15 Z**3 16 T*G**2 17 T*Z**2 18 10**T 19
10**(T*2)
20 G*Z If Teff is the effective temperature, log(G) the decimal
logarithm
of the surface gravity, T and G are defined as: T=log(Teff)-3.7 and
G=log(G)-4.44.
Z is [Fe/H]. Each wavelength point is fittedt to the 20 functions above
and each line in 99010 are the coefficients of these functions. This is
a preliminary model and the final version will be described in the
paper
presenting TGMET.
LL_ELODIE/99011
Internal TGM model Same format as LL_ELODIE/99009 but the coefficient
are
fitted on the whole archive using the atmospheric parameters obtained
from
a fit to LL_ELODIE/99010.