MIDAS'da Sıralı Tayf Analizi
Uygulaması
- Hasan H. Esenoğlu
- İ.Ü. Fen Fakültesi Astronomi ve Uzay Bilimleri Bölümü
Tayf Dalgaboyu İndirgemesi
inmidas
- Midas0001>intape/fits 1-39 bec
tape1
- Midas0002>creat/disp ? 400,600,1,1
- Midas0003>load/lut rainbow
- Midas0004>load/ima
bec0029
.......................................................................................
(yıldız)
- Midas0005>get/curs ................................
(4
nokta tayf kenarlarından, 2 de uzaydan değer alınır)
- Midas0006>average/col
tr051 = bec0029 @0186,@0195
........................................... (yıldız)
- Midas0007>average/col tr052 =
bec0029 @0201,@0237 ............................................. (uzay)
- Midas0008>creat/graph ? 900,350,1,1
- Midas0009>plot/row tr052
- Midas0010>comput/ima str051
= (tr051-tr052) * 9 ...................................................... (yıldız)
- Midas0011>comput/ima str052 =
tr052 * 9 ...................... ( tayfın genişliği:195-186=<9>;
uzay)
- Midas0012>write/desc str051
ident "QU Vul 6931TR"
- Midas0013>write/desc str052 ident
"QU Vul 6931SK"
- Midas0014>plot/row str051
- Midas0015>load/ima bec0027
............................................................................
(standart
yıldız)
- Midas0016>get/curs
............................... (4
nokta tayfın kenarlarından, 2 de uzaydan değer alınır)
- Midas0017>average/col
tr061 = bec0027 @0186,@0195
................................ (standart yıldız)
- Midas0018>average/col tr062 =
bec0029 @0149,@0174 ............................................. (uzay)
- Midas0019>plot/row tr062
- Midas0020>modify/gcurs tr062 1
400,600 .............. (CCD'den kaynaklanan gürültü varsa kesilir)
- Midas0021>plot/row tr062
- Midas0022>comput/ima str063 =
(tr061-tr062) * 9 .......................................... (standart yıldız)
- Midas0023>comput/ima str064 =
tr062 * 9 ....................... (tayfın genişliği:195-186=<9>;
uzay)
- Midas0024>load/ima
bec0028
......................................................................................
(Fe-Ar)
- Midas0025>get/curs .................................................................................
(2
nokta tayftan alınır)
- Midas0026>average/col
str065 = bec0028 @0067,@0252
.......................................... (Fe-Ar)
- Midas0027>write/desc
str063 ident "BD+40 40 32 692931TR"
- Midas0028>write/desc str064 ident
"BD+40 40 32 692931SK"
- Midas0029>write/desc str065 ident
"Fe-Ar 692930 CONF"
- Midas0030>plot/row str065
- Midas0031>cut/ima str065 0,4000
- Midas0032>modify/gcurs str065 0,400
................................. (burada örneğin, dalgaboyu
belirlenmiş 4200.0, 4769.89 çift pikli çizgilerin ufak bileşeni ve
4545.08-5657.94 ile 5606.73-6032.12 arasındaki belirlenmemiş kuvvetli iki
çizgi kesilir)
- Midas0033>plot/row str065
-
- UNIX/LINUX>creat iron1.dat
- 4200.0
- 4545.08
- 4657.94
- 4764.89
- 4879.90
- 5270.36
- 5328.53
- 5558.7
- 5606.7
- 6032.12
- 6114.92
- 6172.28
- 6416.32
- 6677.28
- 6752.83
- 6871.29
- 7147.04
- 7272.93
- exit
- type iron1.dat. ....................................... (18 adet
labaratuvar Fe Ar tayfının dalgaboyu değerleri)
- Midas0034>comput/ima conf072 =
str065+0 ...............................................................
(Fe-Ar)
- Midas0035>plot/row conf072
- Midas0036>cut/ima conf072 0,2000
- Midas0037>plot/row conf072
- Midas0038>cut/ima conf072 500,3000
- Midas0039>modify/gcurs conf072 1
1,400
- Midas0040>create/tab iron1 1 18
iron1.dat. ......................................... (iron1.tb;1'i oluşturur)
- Midas0041>name/col iron1.dat #1 :wave
- Midas0042>read/tab iron1
- Midas0043>set/context
spec ................................................................. (veya
set/context long)
- Midas0044>search/line
conf072 5,100 ? gauss
- Midas0045>set/graph pmode=1
- Midas0046>plot/row conf072
- Midas0047>iden/line
.............................................. (veya iden/long conf001 ystart=0 iron1 4;
- .................................................................................................
buna
göre iron1.tb;1'i oluşturur)
- iden ? :4200.0 ¿
- iden ? :4545.08 ¿
- iden ? :4764.89 ¿
- iden ? :4879.90 ¿
- iden ? :5558.7 ¿
- iden ? :6032.12 ¿
- iden ? :6416.32 ¿
- iden ? :6677.28 ¿
- iden ? :7147.04 ¿
- iden ? :7272.93 ¿
- (Bir tane boş yer işaretlenir)
- Midas0048>calib/line 2,4 ? iron1
................................... (veya calib/long 2 4 iron1.tb;1
gauss)
- Midas0049>get/gcurs
.............................................. (Eğer
çok değişmiş noktalar varsa bunlar
- ...........................................................................
işaretlenerek
hangi çizgiler olduğu tespit edilir)
- Midas0050>creat/ima
ref082 1,1000 3500,4
- Midas0051>rebin/wave str065 fear92
ref082 ................................................................(Fe-Ar)
- Midas0052>plot/row fear92
- Midas0053>cut/ima fear92 100,2000
- Midas0054>plot/row fear92
- Midas0055>center/gauss gcurs ? em
.................................... (Çizgilerin
iki kanadı işaretlenerek
- ...........................................................................
dalgaboyu
kalibrasyonunun doğruluğu test edilir)
- Midas0056>plot/row
str051
- Midas0057>rebin/wave str051 str066
ref082
.................................................................. (yıldız)
- Midas0058>write/desc str066 ident
"QU VUL 6931TRWL"
- Midas0059>plot/row str066
- Midas0060>center/gauss gcurs ? em ..................... (Yıldızın
tayfında belirlenmemiş çizgilerin iki
- ......................................................kanadı işaretlenerek gözlenen
dalgaboyu ve FWHM elde edilir)
- Midas0061>set/graph
pmode=2
- Midas0062>plot/row str066
- Midas0063>copy/graph ps
- Midas0064>rebin/wave str052 str067
ref082 .................................................................
(uzay)
- Midas0065>write/desc str067 ident
"QU Vul 6931SKWL"
- Midas0066>plot/row str067
- Midas0067>plot/row str063
- Midas0068>rebin/wave str063 str068
ref068 ................................................. (standart yıldız)
- Midas0069>write/desc str068 ident
"BD+40 40 32 6929 TRWL"
- Midas0070>plot/row str068
- Midas0071>copy/graph ps
- Midas0072>rebin/wave str065 str069
ref082 ............................................................... (Fe-Ar)
- Midas0073>write/desc str069 ident
"Fe-Ar 6930 conf WL"
- Midas0074>plot/row str069
- Midas0075>creat/icat icat str0*.bdf
- Midas0076>read/icat icat
- Midas0077>outtape/fits
icat,14-15,19-20,25-26,30 tape0 an 6250,1 o
- 1 STR051.BDF
QU VUL 16931TR
- 2 STR052.BDF
QU VUL 16931SK
- 3 STR066.BDF
QU VUL
16931TRWL
- 4 STR067.BDF
QU VUL
16931SKWL
- 5 STR063.BDF
BD+40 40 32 16929TR
- 6 STR064.BDF
BD+40 40 32 16929SK
- 7 STR068.BDF
BD+40 40 32 16929TRWL
- [7 adet dosya manyetik teybe veya DAT kasete (okuyucu unitesinin
tanimlanan ismi yazilarak) alınmış olur]
Tayf Akı İndirgemesi
- Midas0001>creat/graph ? 900,350,1,1
- Midas0002>set/graph pmode=1
- Midas0003>plot/row str066
- Midas0004>plot/row str068
- Midas0005>creat/icat icat str0*.bdf
- Midas0006>read/icat icat
-
- UNIX/LINUX>creat bd4040.dat
- 3636. 35.6 49.0
- 3704. 34.6 49.0
- 3862. 40.5 49.0
- 4036. 38.5 49.0
- 4167. 37.1 49.0
- 4255. 35.9 49.0
- 4666. 33.2 49.0
- 4566. 32.3 49.0
- 4785. 29.4 49.0
- 5000. 27.0 49.0
- 5263. 25.5 98.0
- 5556. 22.5 98.0
- 5840. 20.0 98.0
- 6056. 19.4 98.0
- 6436. 17.0 98.0
- 6790. 15.2 98.0
- 7100. 14.3 98.0
- 7550. 11.6 98.0
- exit
- type bd4040.dat. ........................... (18 adet standart yıldızın
dalgaboyuna karşılık akı değerleri)
-
- Midas0007>creat/tab bd4040 3 18
bd4040.dat. .............................. (bd4040.tb;1'i oluşturur)
- Midas0008>name/col bd4040.dat #1 :x
wavelength
- Midas0009>name/col bd4040.dat #2
:flux
- Midas0010>name/col bd4040.dat #3
:bin width
- Midas0011>read/tab bd4040
- Midas0012>set/context
spec ................................................................. (veya
set/context long)
- Midas0013>resp/spec
resp01 str068 bd4040 start n poly 5
- Midas0014>comput/ima str012 =
str068/resp01 ............................................. (standart yıldız)
- Midas0015>write/desc str012 ident
"BD+404032 6929FLUX"
- Midas0016>read/icat icat
- Midas0017>comput/ima star01 =
str066/3.0 ............................................ (yıldızın poz süresi
- .............. standart yıldızın poz süresinin üç katı olduğundan
üçe bolünerek poz süreleri eşitleniyor)
- Midas0018>comput/ima str013
= star01/resp01 ...........................................................
(yıldız)
- Midas0019>write/desc str013 ident
"QU VUL 6931FLUX"
- Midas0020>plot/row str012
- Midas0021>plot/row str013
- Midas0022>set/graph pmode=2
- Midas0023>plot/row str013
- Midas0024>copy/graph ps
- Midas0025>plot/row str012
- Midas0026>copy/graph ps
- Midas0027>read/icat icat
- Midas0028>outtape/fits icat, 12-13
tape0 an 6250,1 o ................ [2 adet dosya manyetik teybe
- ................................. veya DAT kasete (okuyucu unitesinin
tanimlanan ismi yazilarak) alınmış olur]
- Midas0029>plot/axes
3800,5200,500,100 0,4,1,.5 ? "Wavelength" "Flux"
- Midas0030>overplot/row str013
- Midas0031>label/graph "QU Vul
16931 Flux" curs 0 1 1
- Midas0032>copy/graph ps