Kapelari
From MediaWiki
m (→Hindu bias) |
m (→Hindu bias) |
||
Line 348: | Line 348: | ||
|[[File:PGC32768_hindu_mc1.gif]]||[[File:PGC32768_hindu_mc2.gif]]||[[File:PGC32768_hindu_mc3.gif]] | |[[File:PGC32768_hindu_mc1.gif]]||[[File:PGC32768_hindu_mc2.gif]]||[[File:PGC32768_hindu_mc3.gif]] | ||
|} | |} | ||
- | + | ---- | |
- | + | ||
{| class="wikitable" align="center" | {| class="wikitable" align="center" | ||
! snr = 20 bias = 0.1 !! snr = 20 bias = 0.3 !! snr = 20 bias = 0.5 | ! snr = 20 bias = 0.1 !! snr = 20 bias = 0.3 !! snr = 20 bias = 0.5 |
Revision as of 07:19, 1 July 2013
Contents |
SDSS spektri
- Pokrivenost (3800,9200) Å
- Precnik fajbera 3" (180νm)
- Rezolucija = 1500 @ 3800 Å, 2500 @ 9000 Å
- Pixel spacing log-wavelength (10-4dex)
- Jedinice merenja fluksa 10-17 erg cm-2 s-1 Å-1
Struktura SDSS DR9 spektara:
FLUX [FLOAT] LOGLAM [FLOAT] IVAR [FLOAT] AND_MASK [LONG] OR_MASK [LONG] WDISP [FLOAT] SKY [FLOAT] MODEL [FLOAT]
Emisione linije
; [HeI]* [NeV] [NeVI] [OII] [OII] [NeIII]* HeI [NeIII]* He* [SII] Hd Hg lines=[3204.08, 3346.79, 3426.85, 3727.092, 3729.875, 3869.79, 3889.0, 3968.52, 3971.19, 4072.3, 4102.89, 4341.68, $ ; [OIII] HeII* [ArIV]* [ArIV]* Hb [OIII] [OIII] [OIII] [NI]* [NI]* 4364.436, 4687.05, 4712.62, 4741.43, 4862.68, 4932.603, 4960.295, 5008.240, 5199.35, 5201.84, $ ; sky HeI* NaI* sky NaI* sky [OI]* [OI] sky* [OI] [NI] [NII] 5578.5, 5877.23, 5891.63, 5894.6, 5897.63, 6301.7, 6301.94, 6302.046, 6364.76, 6365.536, 6529.03, 6549.86, $ ; Ha [NII] [SII] [SII] [ArV]** sky [ArIII]** [OII]** [OII]** [FeXI]** 6564.61, 6585.27, 6718.29, 6732.67, 7007.6323, 7246.0, 7137.3670, 7321.9167, 7332.2198, 7894.0711, $ ; [SIII]*** [SIII]*** 8832.33, 9071.39]
Linije su date u vakuumu. Konvertovane su iz 'vazduha' koristeci idl-ovu pro airtovac.pro. IAU standardna konverzija data je u Morton (1991, ApJS, 77, 119):
AIR = VAC / (1.0 + 2.735182E-4 + 131.4182 / VAC^2 + 2.76249E8 / VAC^4)
Linije oznacene sa
- * su emisione linje iz Gandalf_V5.1
- ** iz Appenzeller and Ostreicher 1988, AJ 95, 45
- *** iz Osterbrock-a (Table 3.12) (znacajne za sf galaksije)
Druge linije preuzete su iz liste linija SDSS DR7.
Napravila sam proceduru za uklanjanje emisionih linija iz spektara emission_lines.pro, koji poziva fitcont.pro za fitovanje kontinuuma. Program izlistava emsione linije i ukoliko su znacajno iznad kontinuuma i imaju dovoljno tacaka (>4) fituje ih gausijanom i oduzima od ulaznog fluksa. Kod se koristi tako sto se prvo ucita deredsiftovani spektar galaksije koji iz SDSS DR9 fit fajla pravi sdss2txt.pro. Koraci su sledeci:
sdss2txt, 'NGC2410.fits', redshift=0.01546652 ;ciji je izlaz NGC2410_restframe.txt i nalazi se u istom dir gde i fits gal='NGC2410_restframe.txt' & readcol, gal, F='(F,F,F)', wave, flux,fluxivar ; ucitavam spektar out=emission_lines(wave,flux,fluxivar) ; out -- spektar od kojeg su oduzete emisione linije iz liste gore plot,wave,flux,background='FFFFFF'x,col=255 oplot,wave,out,col='000000'x ;rezultat
Zvezdane biblioteke
Vazdekis
Struktura:
- ppxf
- ppxf_kinematics_example_sdss_vazdekis.pro
- sample_obj_pmf_z_sn.csv
- spectra/Rbi1.30z*
- sdssdr8/spec*fits
- log_rebin.pro
- ppxf.pro
- ppxf_determine_goodpixels.pro
Korisceni templejti (4795 - 5465)AA nalaze se u spectra/ i to su templejti koji dolaze sa ppxf kodom. U pitanju su sinteticki spektri. Korisceni SDSS DR8 spektri nalaze se u sdssdr8/. Dodati fajl sample_obj_pmf_z_sn.csv koristi se da se iz njega izdvoji ime galaksije, njen crveni pomak i snr. Skripta vazdekis_dr8.sh pustena je u petlji:
for i in $(ls sdssdr8/spec*fits); do ./vazdekis_dr8.sh $i 2>>1; done > log
i kreirace dodatne fajlove:
- ppxf/vazdekis/veldisp.dat
- ppxf/vazdekis/Figures/*.ps # ps slike fitova
head -1 ppxf/vazdekis/veldisp.dat # galaxy V Vdisp h3 h4 chi2dof V_err Vdisp_err h3_err h4_err bias template NGC3042 58.3 179.9 -0.0310 0.0329 4.e-07 3.3 4.5 0.0169 0.0178 0.7 spectra/Rbi1.30z+0.20t17.78.fits
Rezultati su sacuvani u vazdekis.tgz:
tar tvf vazdekis.tgz vazdekis_dr8.sh ppxf_kinematics_example_sdss_vazdekis.pro ppxf_mcbiaspar.pro sample_obj_pmf_z_sn.csv veldisp.dat log
- σ=160.2 h3=0.001 h4=-0.033 χ2=2e-5 bias=0.7 snr=79
- σ=213.5 h3=-0.009 h4=0.039 χ2=6e-6 bias=0.7 snr=71
- σ=217.5 h3=-0.012 h4=-0.003 χ2=1e-1 bias=0.7 snr=9
Elodie
Spectralni opseg (AA):
- SDSS: 3808 - 9266
- Elodie: 3900 - 6800
Svih 1386 zvezda imaju rupe. Od toga samo 3 imaju rupe vece od 100AA, ali jos 388 imaju rupe vece od 2 piksela sa jednog ili oba kraja, sto ostavlja konacno 998 zvezda sa rupama manjim od 100AA koje nisu na krajevima spektara ( elodie_good_files_ui.dat). One su propustene kroz midas da im ispegla rupe ( midas_elodie.prg koji poziva script_elodie.sh). Finalno, koristi se 998 splajnovanih templejta dobijenih na ovaj nacin.
Glavna skripta elodie_dr8.sh poziva ppxf_kinematics_example_sdss_elodie.pro koristeci sample_obj_pmf_z_sn.csv.
- σ=157.6 h3=0 h4=-0.03 χ2=2e-5 bias=0.4 snr=79
- σ=209.2 h3=0.001 h4=0.001 χ2=6e-6 bias=0.4 snr=71
- σ=270.9 h3=-0.035 h4=0.02 χ2=8e-2 bias=0.4 snr=9
INDO-U.S.
Indo-U.S. Coude Feed Spectral Library: 1273 stars covering the region (3460, 9464) Å at a spectral resolution of 1.35 Å (FWHM), σ ∼ 30 km/s, R ∼ 4200. Treba skinuti tri fajla:
- index.html (sa wget)
- table2.txt
- table3.txt
========== index.html ========== [TXT] 249.fits 10-Feb-2004 15:35 802K [TXT] BD+01_2916.fits 10-Feb-2004 15:44 667K [TXT] BD-01_2582.fits 10-Feb-2004 15:44 397K [TXT] G_4-44.fits 10-Feb-2004 15:46 532K [TXT] NGP_29_129.fits 10-Feb-2004 15:47 532K ========== table2.dat ========== stellar_id coverage gaps 1-10(A10) 51-63(A27) 65-91(A13) ----------------------------------------------------------------------------- G 4-44 02 51 58.36 +11 22 11.9 8.38 0.54 3 3465.0-7265.8 G 12-21 12 12 01.37 +13 15 40.6 10.18 0.47 4 3465.0-8104.2 5413.0-5465.4 BD-11 4126 16 23 26.91 -11 35 40.6 10.40 0.80 4 3465.0-8107.4 5415.4-5462.6 BD+01 2916 14 21 45.26 +00 46 59.2 9.65 1.28 4 3465.0-8647.0 HD 249 00 07 22.56 +26 27 02.2 7.33 1.00 5 3465.0-9469.0 6789.4-6810.6 NGP 29 129 12 39 47.00 +28 19 00.0 11.13 0.85 3 3465.0-7312.6 ----------------------------------------------------------------------------- ========== table3.dat ========== stellar_id spec type 1-10(A10) 12-21(A10) ----------------------------------------------------------------------------------- G 4-44 G5 G5V 6.2 5750 4.11 -0.69 2000A&A...353..722N (Nissen) BD-11 4126 K3V K3V 5000 4.30 0.20 1997AJ....114..376C (Castro) BD+09 3223 Flat 67.0 5350 2.00 -2.26 2000ApJ...544..302B (Burris) HD 358 B8IVmnp... B8V -11.7 13622 4.00 -0.30 1979ApJS...41..675H (Heacox) NGP 29 129 K1V K2V 4.00 0.30 1985AJ.....90..803R (Rose) -----------------------------------------------------------------------------------
1. Check if gaps from the headers coincide with the gaps from the table2.txt file: There are 721 stars WITHOUT GAPS, i.e. 1273-721=552 with gaps.
cut -c65-91 table2.dat | sed 's/^ * $/null/g' | grep -w null | wc -l
- I had to download the whole lib to test headers of the fits files to confirm that there are 721 stars WITHOUT GAPS, i.e. 552 with gaps accord to the headers.
for i in $(ls *fits); do gaps=$(gethead gaps $i); if ! [ -z "$gaps" ]; then echo "$i: $gaps"; fi;done | wc -l
2. Make wget script using only those stars with no gaps and good wave range (excuding stars with gaps and insufficient wave range):
# SPECTRAL TYPE: cut -c11-21 table3.dat | tr -s [A-Za-z] | tr -d : | tr -d . | tr -d + | tr -d / | tr '-' '_' > spTypes.txt # URLS: awk -v baseurl="http://www.noao.edu/cflib/V1/BINTABLE/" '/fits/{print baseurl$2}' index.html > all_spectra.txt # NAMES: cut -c1-10 table2.dat | sed s/HD//g | awk '{OFS="_";print $1,$2,$3}' | sed 's/__$//g' | sed 's/_$//g' > names.dat paste names.dat coverage.dat gaps_null.dat spTypes.txt | grep null | awk '($2<3808 && $3>9266){print}' | wc -l
FINALY, there are 441 stars in the SDSS wave range [3808:9266] without gaps.
3. Rename downloaded stars according to their spec type and stellar identification
paste names.dat coverage.dat gaps_null.dat spTypes.txt | grep null | awk '($2<3808 && $3>9266){print}' | \ awk -v baseurl="http://www.noao.edu/cflib/V1/BINTABLE/" '{print "wget -O hindu_"$5"_"$1".fits", baseurl$1".fits"}' \ > hindu_download_NEW.sh # Example: wget -O hindu_G3V_G_12-24.fits http://www.noao.edu/cflib/V1/BINTABLE/G_12-24.fits
4. Determine spec coverage of those - make it the same and widest as possible (midas_hindu.prg) i.e. (3780,9270)AA
for i in $(ls hindu*[0-9].fits); do ./script_hindu.sh $i; done
- σ=167.8 h3=0.008 h4=-0.018 χ2=2e-5 bias=0.3 snr=79
- σ=328.6 h3=-0.001 h4=0.069 χ2=8e-6 bias=0.3 snr=71
- σ=327.2 h3=0.000 h4=0.008 χ2=1e-1 bias=0.3 snr=9
Ako se hindu templejti seku kao Vazdekis [4795:5465], dobijaju se dosta bolji rezultati:
galaxy | σ | h3 | h4 | χ2 |
---|---|---|---|---|
NGC2410 | 158.3 | -0.003 | -0.074 | 0.000 |
UGC4956 | 201.7 | -0.004 | 0.003 | 0.000 |
NGC2672 | 204.9 | -0.017 | -0.003 | 0.098 |
Stelib
- The current version of the STELIB library contains the spectra of the 255 stars
- The spectra are corrected for interstellar extinction and radial velocity.
- The initial dispersion has been modified to 1A/pixel.
- The files contain 6700 pixels from 3200AA to 9900AA.
Sa Vizijera iz ReadMe fajla nasla sam table.dat sa informacijom o spektralnom tipu zvezda.
- 1. Preimenovati skinute fits fajlove prema nj. spektralnom tipu:
cut -c1-11 -c36-46 -c124-143 table.dat | tr -d / | tr -d =:= | tr -d =+= | tr -d . | \ sed -e 's/fits/.fits/g' -e 's/V32/V3.2/g' | awk '! ($4 ~ /^$/){print "ln -s "$4,"stelib_"$3"_"$1$2".fits"}' \ > make_links.sh # ex. stelib_A2V_HD097633.fits -> HD097633_V3.2.fits # There is one missing : LINE=68 in table.dat -- no download link # SO THERE ARE 254 stellar templates
- 2. Naci opseg talasnih duzina, da bi bile uzete samo one zvezde koje obuhvataju SDSS opseg [3808:9206]AA.
for i in $(ls stelib*fits); do echo $i $(./range.sh $i); done > ranges.dat awk '{print $2}' ranges.dat | sort -n | tail -1 max of the lower limit: 3204.3520 awk '{print $3}' ranges.dat | sort -nr | tail -1 min of the upper limit: 9887.0505 # SO THE RANGE IN COMMON IS [3205:9887] # To explore the full SDSS range templates should be cut in [3780:9270]AA.
- 3. Ima zvezda sa 'rupama' (flux=0) koje uzrokuju ppxf da pukne sa porukom too WIERD distribution. Odlucila sam da odbacim samo one koje imaju vise od jednog defektog piksela na krajevima i vecu rupu od 100AA. Ovo ispituje procedura checkNans.pro. Ona pravi stelib_gaps_edges.dat koji lista ime zvezde, ima li rupa (yes/no), najvecu rupu u pikselima i angstremima i, konacno medijanu (za kasniju normalizaciju spektra u midas-u).
Primer: stelib_sdG2_HD064090.fits yes 446 pixels 446.3 angstroms 1.31649e-12
# Spektri koji nemaju rupa na krajevima, a izmedju manjih od 100AA ima 109!!! awk '($3==1 && $6<100){print }' stelib_gaps_edges.dat | wc -l # Od njih pravim listu tzv. dobrih fajlova za ppxf. awk '($3==1 && $6<100){print $1}' stelib_gaps_edges.dat > stelib_good_files.dat # Ubacujem ih u skriptu koja zapravo samo poziva midas da ih isece i splajnom ispegla rupe (midas_stelib.prg) while read line; do ./script_stelib.sh $line ; done < stelib_good_files.dat
- σ=133.3 h3=0.028 h4=-0.197 χ2=7e-5 bias=0.3 snr=79
- σ=153.3 h3=0.000 h4=-0.001 χ2=4e-5 bias=0.3 snr=71
- σ=209.9 h3=0.001 h4=-0.001 χ2=1e-1 bias=0.3 snr=9
Ako se stelib templejti seku kao Vazdekis [4795:5465], dobijaju se dosta bolji rezultati:
galaxy | σ | h3 | h4 | χ2 |
---|---|---|---|---|
NGC2410 | 146.4 | 0.011 | -0.078 | 0.000 |
UGC4956 | 186.5 | -0.014 | 0.010 | 0.000 |
NGC2672 | 213.4 | -0.016 | 0.003 | 0.099 |
Miles
- Biblioteka sadrzi 985 spektara
- Spektri snimljeni 2.5m INT teleskopom
- Spektraln pokrivenost [3525-7500]Å
- Spektralna rezolucija 2.50Å (FWHM)
- Trenutna verzija MILES_library_v9.1.tar.gz
Sa Vizier-a sam skinula catalog.dat sa nazivima fajlova i tipovima zvezda.
cut -c59-72 catalog.dat | sed 's/^ * $/null/g' > type.dat cut -c1-5 catalog.dat | sed -e s/s/m/g -e s/$/V/g > names.dat
Pomocu checkNans.pro napravila sam txt fajlove od onih zvezda koje na krajevima nemaju rupe vece od 1AA i cija najveca rupa nije veca od 100AA, popunila rupe Nan-ovima (da bi ih midas ispeglao) i isekla ih da pokriju SDSS opseg [3800.6:7409.6]. Ovako pravim 979 fajlova, koje ce dalje peglati midas.
paste names.dat type.dat | awk '{print "if [ -e "$1".txt ]; then ln -s",$1".txt","miles_"$1"_"$2".txt; fi"}' > tmp cat tmp | tr -d : | tr -d + | tr -d ? | tr -d "..." | tr -d / | tr -d \( | tr -d \) | sed s/txt/.txt/g | sed s/-/_/g | \ sed s/_e/-e/g | sed s/_s/-s/g > make_links.sh head -1 make_links.sh # if [ -e m0001V.txt ]; then ln -s m0001V.txt miles_m0001V_G3.txt; fi
Dakle, nekoliko ih treba ispeglati (popuniti rupe splajn interpolacijom), a ostale prosto pretvoriti u fits. Ovo radi midas_miles.pro:
for i in $(ls miles_m0*.txt); do sed s/infile/${i%.*}/ midas_miles.prg > midas_${i%.*}.prg; inmidas < midas_${i%.*}.prg; \ rm midas_${i%.*}.prg; done
Konacno, ima 979 fajlova naziva miles_m0001V_G3_spline.fits, koje se daju miles_dr8.sh koja poziva ppxf_kinematics_example_sdss_miles.pro. Defoltni bias koji izracunava ppxf bias=0.3.
- σ=161.0 h3=0.0017 h4=-0.022 χ2=2e-5 bias=0.3 snr=79
- σ=217.2 h3=0.020 h4=0.040 χ2=5e-6 bias=0.3 snr=71
- σ=341.3 h3=-0.009 h4=0.012 χ2=1e-1 bias=0.3 snr=9 (Sa bias=0.7 dobija se σ=353.0 h3=-0.004 h4=0.007 χ2=1.e-01)
Test bias parameter
Vazdekis bias
Vazdekis templejti za snr = 10 (UGC10498), snr = 20 (PGC10766) i snr = 30 (PGC58707) za razlicite bajase.
- snr = 10 AND bias = 0.1, 0.2, 0.3
snr = 10 bias = 0.1 | snr = 10 bias = 0.2 | snr = 10 bias = 0.3 |
---|---|---|
snr =20 bias = 0.1 | snr = 20 bias = 0.3 | snr = 20 bias = 0.5 |
snr = 30 bias = 0.1 | snr = 30 bias = 0.3 | snr =30 bias = 0.5 |
Elodie bias
- snr = 10 AND bias = 0.1, 0.2, 0.3
snr = 10 bias = 0.1 | snr = 10 bias = 0.2 | snr = 10 bias = 0.3 |
---|---|---|
snr = 10 bias = 0.1 | snr = 10 bias = 0.2 | snr = 10 bias = 0.3 |
- snr = 20 AND bias = 0.1, 0.3, 0.5
snr = 20 bias = 0.1 | snr = 20 bias = 0.3 | snr = 20 bias = 0.5 |
---|---|---|
- snr = 30 AND bias = 0.1, 0.3, 0.5
snr = 30 bias = 0.1 | snr = 30 bias = 0.3 | snr = 30 bias = 0.5 |
---|---|---|
Hindu bias
- snr = 10 AND bias = 0.1, 0.2, 0.3
snr = 10 bias = 0.1 | snr = 10 bias = 0.2 | snr = 10 bias = 0.3 |
---|---|---|
snr = 10 bias = 0.1 | snr = 10 bias = 0.2 | snr = 10 bias = 0.3 |
snr = 20 bias = 0.1 | snr = 20 bias = 0.3 | snr = 20 bias = 0.5 |
---|---|---|
- snr = 30 AND bias = 0.1, 0.3, 0.5
snr = 30 bias = 0.1 | snr = 30 bias = 0.3 | snr = 30 bias = 0.5 |
---|---|---|
Miles bias
- snr = 10 AND bias = 0.1, 0.2, 0.3
- snr = 20 AND bias = 0.1, 0.3, 0.5
- snr = 30 AND bias = 0.1, 0.3, 0.5
Rezultati
- Grafici crtani uz uslov da su sve disperzije pozitivne i manje od 400. Iskljuceno ukupno 23 galaksije, od cega 14 jer u SDSS bazi nema podataka za disperziju. Ostale zbog nedostatka templejta.
Ako se hindu, stelib i miles templejti iseku kao Vazdekis [4795:5465], dobijaju se dosta bolji rezultati:
galaxy | σ_vaz | σ_hin | σ_ste | σ_mil | h3_vaz | h3_hin | h3_ste | h3_mil | h4_vaz | h4_hin | h4_ste | h4_mil |
---|---|---|---|---|---|---|---|---|---|---|---|---|
NGC2410 | 160.2 | 158.3 | 146.4 | 164.4 | 0.001 | -0.003 | 0.011 | 0.014 | -0.033 | -0.074 | -0.078 | -0.030 |
UGC4956 | 213.5 | 201.7 | 186.5 | 201.1 | -0.009 | -0.004 | -0.014 | 0.005 | 0.039 | 0.003 | 0.010 | 0.005 |
NGC2672 | 217.5 | 204.9 | 213.4 | 204.9 | -0.012 | -0.017 | -0.016 | -0.009 | -0.003 | -0.003 | 0.003 | -0.004 |
Inace, svi su radjeni sa bias=0.7, sto je defoltni bias za Vazekis stelib.
Poredjenje Hidu u celom opsegu (hin_full) i Miles i Hindu templejta secenih u zajednickom (Miles=3800:4700) opsegu sa bias=0.7:
galaxy | σ_hin_full | σ_hin | σ_mil | h3_hin_full | h3_hin | h3_mil | h4_hin_full | h4_hin | h4_mil |
---|---|---|---|---|---|---|---|---|---|
NGC2410 | 167.3 | 158.3 | 160.6 | 0.010 | 0.011 | 0.013 | -0.017 | -0.022 | -0.019 |
UGC4956 | 240.3 | 224.7 | 217.6 | -0.0005 | 0.001 | 0.017 | 0.053 | 0.046 | 0.034 |
NGC2672 | 294.7 | 396.4 | 356.7 | -0.0005 | -0.017 | -0.005 | 0.000 | 0.012 | 0.008 |