Kapelari

From MediaWiki

Revision as of 14:26, 9 July 2013 by Ana (Talk | contribs)
Jump to: navigation, search

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

image:idl.gif

Zvezdane biblioteke

Vazdekis

Struktura:


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
==========
[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
NGC2410158.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
NGC2410146.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 bias = 0.1 snr = 10 bias = 0.2 snr = 10 bias = 0.3
image:UGC10498_mc1.gifimage:UGC10498_mc2.gifimage:UGC10498_mc3.gif
snr =20 bias = 0.1 snr = 20 bias = 0.3 snr = 20 bias = 0.5
File:PGC10766_mc1.gifFile:PGC10766_mc3.gifFile:PGC10766_mc5.gif
snr = 30 bias = 0.1 snr = 30 bias = 0.3 snr =30 bias = 0.5
File:PGC58707_mc1.gifFile:PGC58707_mc3.gifFile:PGC58707_mc5.gif

Elodie bias


snr = 10 bias = 0.1 snr = 10 bias = 0.2 snr = 10 bias = 0.3
File:UGC10498_elodie_mc1.gifFile:UGC10498_elodie_mc2.gifFile:UGC10498_elodie_mc3.gif
snr = 10 bias = 0.1 snr = 10 bias = 0.2 snr = 10 bias = 0.3
File:PGC32768_elodie_mc1.gifFile:PGC32768_elodie_mc2.gifFile:PGC32768_elodie_mc3.gif
snr = 20 bias = 0.1 snr = 20 bias = 0.2 snr = 20 bias = 0.3
File:PGC10766_elodie_mc1.gifFile:PGC10766_elodie_mc2.gifFile:PGC10766_elodie_mc3.gif
snr = 20 bias = 0.4 snr = 20 bias = 0.5 snr = 20 bias = 0.6
File:PGC10766_elodie_mc4.gifFile:PGC10766_elodie_mc5.gifFile:PGC10766_elodie_mc6.gif
snr = 30 bias = 0.1 snr = 30 bias = 0.2 snr = 30 bias = 0.3
File:PGC58707_elodie_mc1.gifFile:PGC58707_elodie_mc2.gifFile:PGC58707_elodie_mc3.gif
snr = 30 bias = 0.4 snr = 30 bias = 0.5 snr = 30 bias = 0.6
File:PGC58707_elodie_mc4.gifFile:PGC58707_elodie_mc5.gifFile:PGC58707_elodie_mc6.gif
snr = 40 bias = 0.1 snr = 40 bias = 0.3 snr = 40 bias = 0.5
File:NGC223_elodie_mc1.gifFile:NGC223_elodie_mc3.gifFile:NGC223_elodie_mc5.gif
snr = 40 bias = 0.6 snr = 40 bias = 0.7 snr = 40 bias = 0.9
File:NGC223_elodie_mc6.gifFile:NGC223_elodie_mc7.gifFile:NGC223_elodie_mc9.gif
snr = 50 bias = 0.3 snr = 50 bias = 0.5 snr = 50 bias = 0.6
File:NGC5831_elodie_mc3.gifFile:NGC5831_elodie_mc5.gifFile:NGC5831_elodie_mc6.gif
snr = 50 bias = 0.7 snr = 50 bias = 0.8 snr = 50 bias = 0.9
File:NGC5831_elodie_mc7.gifFile:NGC5831_elodie_mc8.gifFile:NGC5831_elodie_mc9.gif
snr = 60 bias = 0.1 snr = 60 bias = 0.3 snr = 60 bias = 0.5
File:NGC2424_elodie_mc1.gifFile:NGC2424_elodie_mc3.gifFile:NGC2424_elodie_mc5.gif
snr = 60 bias = 0.7 snr = 60 bias = 0.9 snr = 60 bias = 1.0
File:NGC2424_elodie_mc7.gifFile:NGC2424_elodie_mc9.gifFile:NGC2424_elodie_mc10.gif
snr = 70 bias = 0.3 snr = 70 bias = 0.5 snr = 70 bias = 0.6
File:NGC4161_elodie_mc3.gifFile:NGC4161_elodie_mc5.gifFile:NGC4161_elodie_mc6.gif
snr = 70 bias = 0.7 snr = 70 bias = 0.8 snr = 70 bias = 0.9
File:NGC4161_elodie_mc7.gifFile:NGC4161_elodie_mc8.gifFile:NGC4161_elodie_mc9.gif

Hindu bias

snr = 10 bias = 0.1 snr = 10 bias = 0.2 snr = 10 bias = 0.3
File:UGC10498_hindu_mc1.gifFile:UGC10498_hindu_mc2.gifFile:UGC10498_hindu_mc3.gif
snr = 10 bias = 0.1 snr = 10 bias = 0.2 snr = 10 bias = 0.3
File:PGC32768_hindu_mc1.gifFile:PGC32768_hindu_mc2.gifFile:PGC32768_hindu_mc3.gif
snr = 20 bias = 0.1 snr = 20 bias = 0.2 snr = 20 bias = 0.3
File:PGC10766_hindu_mc1.gifFile:PGC10766_hindu_mc2.gifFile:PGC10766_hindu_mc3.gif
snr = 20 bias = 0.4 snr = 20 bias = 0.5 snr = 20 bias = 0.6
File:PGC10766_hindu_mc4.gifFile:PGC10766_hindu_mc5.gifFile:PGC10766_hindu_mc6.gif
snr = 30 bias = 0.1 snr = 30 bias = 0.3 snr = 30 bias = 0.5
File:NGC2424_hindu_mc1.gifFile:NGC2424_hindu_mc3.gifFile:NGC2424_hindu_mc5.gif
snr = 40 bias = 0.3 snr = 40 bias = 0.4 snr = 40 bias = 0.5
File:NGC223_hindu_mc3.gifFile:NGC223_hindu_mc4.gifFile:NGC223_hindu_mc5.gif
snr = 40 bias = 0.6 snr = 40 bias = 0.7 snr = 40 bias = 0.9
File:NGC223_hindu_mc6.gifFile:NGC223_hindu_mc7.gifFile:NGC223_hindu_mc9.gif
snr = 50 bias = 0.3 snr = 50 bias = 0.5 snr = 50 bias = 0.6
File:NGC5831_hindu_mc3.gifFile:NGC5831_hindu_mc5.gifFile:NGC5831_hindu_mc6.gif
snr = 50 bias = 0.7 snr = 50 bias = 0.8 snr = 50 bias = 0.9
File:NGC5831_hindu_mc7.gifFile:NGC5831_hindu_mc8.gifFile:NGC5831_hindu_mc9.gif
snr = 60 bias = 0.1 snr = 60 bias = 0.3 snr = 60 bias = 0.5
File:NGC2424_hindu_mc7.gifFile:NGC2424_hindu_mc9.gifFile:NGC2424_hindu_mc10.gif
snr = 60 bias = 0.7 snr = 60 bias = 0.9 snr = 60 bias = 1.0
File:PGC58707_hindu_mc1.gifFile:PGC58707_hindu_mc3.gifFile:PGC58707_hindu_mc5.gif
snr = 70 bias = 0.3 snr = 70 bias = 0.5 snr = 70 bias = 0.6
File:NGC4161_hindu_mc3.gifFile:NGC4161_hindu_mc5.gifFile:NGC4161_hindu_mc6.gif
snr = 70 bias = 0.7 snr = 70 bias = 0.8 snr = 70 bias = 0.9
File:NGC4161_hindu_mc7.gifFile:NGC4161_hindu_mc8.gifFile:NGC4161_hindu_mc9.gif

Miles bias

snr = 10 bias = 0.1 snr = 10 bias = 0.2 snr = 10 bias = 0.3
File:UGC10498_miles_mc1.gifFile:UGC10498_miles_mc2.gifFile:UGC10498_miles_mc3.gif
snr = 10 bias = 0.1 snr = 10 bias = 0.2 snr = 10 bias = 0.3
File:PGC32768_miles_mc1.gifFile:PGC32768_miles_mc2.gifFile:PGC32768_miles_mc3.gif
snr = 20 bias = 0.1 snr = 20 bias = 0.2 snr = 20 bias = 0.3
File:PGC10766_miles_mc1.gifFile:PGC10766_miles_mc2.gifFile:PGC10766_miles_mc3.gif
snr = 20 bias = 0.4 snr = 20 bias = 0.5 snr = 20 bias = 0.6
File:PGC10766_miles_mc4.gifFile:PGC10766_miles_mc5.gifFile:PGC10766_miles_mc6.gif
snr = 30 bias = 0.1 snr = 30 bias = 0.2 snr = 30 bias = 0.3
File:PGC58707_miles_mc1.gifFile:PGC58707_miles_mc2.gifFile:PGC58707_miles_mc3.gif
snr = 30 bias = 0.4 snr = 30 bias = 0.5 snr = 30 bias = 0.6
File:PGC58707_miles_mc4.gifFile:PGC58707_miles_mc5.gifFile:PGC58707_miles_mc6.gif
snr = 40 bias = 0.2 snr = 40 bias = 0.3 snr = 40 bias = 0.4
File:NGC223_miles_mc2.gifFile:NGC223_miles_mc3.gifFile:NGC223_miles_mc4.gif
snr = 40 bias = 0.5 snr = 40 bias = 0.6 snr = 40 bias = 0.7
File:NGC223_miles_mc5.gifFile:NGC223_miles_mc6.gifFile:NGC223_miles_mc7.gif
snr = 50 bias = 0.3 snr = 50 bias = 0.5 snr = 50 bias = 0.6
File:NGC5831_miles_mc3.gifFile:NGC5831_miles_mc5.gifFile:NGC5831_miles_mc6.gif
snr = 50 bias = 0.7 snr = 50 bias = 0.8 snr = 50 bias = 0.9
File:NGC5831_miles_mc7.gifFile:NGC5831_miles_mc8.gifFile:NGC5831_miles_mc9.gif
snr = 60 bias = 0.1 snr = 60 bias = 0.3 snr = 60 bias = 0.5
File:NGC2424_miles_mc1.gifFile:NGC2424_miles_mc3.gifFile:NGC2424_miles_mc5.gif
snr = 60 bias = 0.7 snr = 60 bias = 0.9 snr = 60 bias = 1.0
File:NGC2424_miles_mc7.gifFile:NGC2424_miles_mc9.gifFile:NGC2424 miles mc10.gif
snr = 70 bias = 0.3 snr = 70 bias = 0.5 snr = 70 bias = 0.6
File:NGC4161_miles_mc3.gifFile:NGC4161_miles_mc5.gifFile:NGC4161_miles_mc6.gif
snr = 70 bias = 0.7 snr = 70 bias = 0.8 snr = 70 bias = 0.9
File:NGC4161_miles_mc7.gifFile:NGC4161_miles_mc8.gifFile:NGC4161_miles_mc9.gif

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.
Disperzije za sve templejte sa koeficijentima korelacije. Crvena linija y=x.
Disperzija za 172 galaksije nadjene u Ledi. Hindu dat sa bias=0.5. Crvena linija y=x. Donji panel - koeficijenti korelacije.
Crvena linija y=x, a plave prolaze kroz nulu. Histogrami dati za sve koriscene zvezdane biblioteke.

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
NGC2410160.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
Personal tools