Page updated the 2012-05-04
Context:
August 2010, at OHP spectral meeting, we have found criteria that may impact on EW's computation and we then defined some conventions about how to compute EW.
Possible criteria that may impact on EW's computation:
EW's definition is:
![]()
I(λ) : spectrum's intensity at wavelength λ
Ic(λ) : continuum's intensity at wavelength λ
EW's computation needs continuum's intensities Ic which can be extracted from spectrum by several ways. So, it seems that there is no need to normalize the spectrum.
We have found the following criteria that should impact on EW's computation:
- Wavelength bandwidth for integration;
- Continuum extraction method;
- Star's behaviour of course;
- For the same date, EW depends on the observer (experimental bias, spectral reduction, ...);
- Depends on software's version (id. algorithm used for computation).
Here are a pictures of white board used at OHP 2010 meeting workshop on EW:
Solutions to fixe criteria:
Workshop dealed with all points adn conclued to:
- Wavelength bandwidth for integration:
We have to fix wavelength's bandwidth for each stars. If unfixed, 6535-6595 is used. This bandwidth has to include all shapes of the line studed including wings. Wings' limits have to be discussed for each stars and depends on each spectrum's resolution (wider for low resolution spectra).- Continuum extraction method:
We have to choose the continuum extraction method and apply it to each spectra. This is the main source of EW's dispersion with intrinsic star's activity variations.
The continuum extracted needs to follow global slop of left edge and right edge stellar's continuum shape and needs to have the smallest curvature.- Star's behaviour of course:
Stars' activity may vary through hours to mounths.- For the same date, EW depends on the observer:
If preprocessing is well done and all EW's computation is made by only one observer, EW won't vary. Only weather conditions could then influence EW value.- Depends on software's version:
Software's version and name have to be present in fits file's header in order to know witch algorithm have been used to compute EW and continuum. Software's EW computation result has to be compare one time with MIDAS's result.It appears that the main reason for EW's variations is the continuum computed used for EW.
Continuum's extraction:
1. How to choose a continuum?
We can see bellow few continuum extracted from the same star and the EW obtained with it. Computation are done with SpcAudace's spc_ew command on this CI Cyg spectrum.
Continuum has a great impact on EW's value. So whitch continuum is the best?
With this star, multi linear pieces smoothed with factor 1000000 seems to be the best continuum fitting which takes care of global slop of left edge and right edge stellar's continuum and that has the smallest curvature.
But with an other star (the Crb spectrum), low pass filter with deg 2 polynom would be the best:Command: spc_ew the_crb_20100315.fit 6535 6595 6 n 2
Low pass filter with deg 2 polynom that matches Ha line wing's edges.
Obtained with command: spc_ew the_crb_20100315.fit 6535 6595 6 n 2Here is how SpcAudace works in automated mode:
- It extracts continuum with non-linear methods (left picture);
- It finds integrations limits by finding intersection between special smoothed spectrum and extracted continuum (right picture).
2. What makes continuum's extraction difficult?
Continuum extraction's main difficulty comes from continuum's SNR. In the last example, SNR was very low and algorithms have few possibilities how to find continuum median position in wide spread distribution of intensities.
Command: spc_ew the_crb_20100315.fit 6535 6595 6 n 2
However, when SNR is good, all algorithms find closed EW's values for the same bandwidth (6535-6595) on this zeta Tau spectrum:
3. Possible work after continnum's extraction method is defined for a star:
So, before beginning before launching a time study of EW's evolution, it is needed to select the more adapted method for continuum extraction witch depends on the star. After what a pipeline may be applyed.
Here is an example of such pipeline for EW time evolution:
EW evolution from 2010/02/05 to 2010/04/04.
Obtained with command: spc_ewcourbe 6535 6595EW's period computation by fiting with unscaled
Levenberg-Marquardt algorithm.A periodogramm (spc_periodogram SpcAudace command) may resolve the period too if there are anougth points:
Conclusions:
Finally, all criteria are fixed but continuum extraction is the critical step:
- Wavelength bandwidth has to be adapted to the star and spectrum's resolution;
- EW computation must be given with an error estimation;
- Continuum extraction must be done by an algorithm in order to avoid manual bias and time waste, and then makes available scripting for time evolution studies with graphs etc.;
- If star's behaviour doesn't change, the same continuum's extraction method will be applied;
- With cautions described below, starting choices can be shared with the whole survey's team and measurements can be done by anyone without any bias. Scientific tasks can then be shared.
With such methodoly, studing stars would be rather efficiency and sure.
BeSS database and dedicated surveys such as novae monitoring will then give a great amount of information.